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ABSTRACT 


The modeling of picosecond pulse propagatiun on microwave integrated circuit 
interconnections is considered. Autoregressive moving-average (ARMA) and autore- 
gressive (AR) parametric models are derived for lossy dispersive microstrip transmis- 
sion lines and cascaded microstrip step discontinuities. We formulated mathematical 
expressions to relate the model parameters to the physical microstrip properties. 

New lumped-distributed equivalent circuit models are presented. Dispersive pulse 
propagation on high-frequency integrated circuit interconnections is modeled using 
frequency-dependent lumped parameters and lossy distributed transmission-line sec- 
tions. We verified the equivalent circuit models through computer simulations and 
experimental measurements. Modern parameter estimation techniques are applied 
to system identification modeling. We develop several algorithms to estimate the 
model parameters from input and/or output measurements. The pertormance of the 


algorithms are evaluated using computer simulations and experimental results. 
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I. INTRODUCTION 


Microwave and high-speed very large scale integrated (WLSI) circuits are lim- 
ited by the propagation characteristics of the on-chip interconnections [Ref. 1]. The 
effects of signal dispersion and loss become even more predominant as devices and 
circuits are scaled to smaller dimensions. At microwave frequencies, the intercon- 
nections between elements on a dielectric substrate, such as silicon (Si) or gallium 
arsenide (GaAs) where considerable wavelength reductions occur, must be treated as 
planar waveguide structures. The analysis and design of circuits consisting of these 
guided wave structures are facilitated by the use of equivalent circuits. Accurate 
high frequency characterization of picosecond pulse propagation on these structures 
including dispersion and losses requires extensive numerical techniques. The em- 
pirical equations that have been derived for the propagation parameters from these 
numerical results do not lead to an equivalent circuit model with realizable series and 
shunt branches. Furthermore, reliable and accurate empirical models for many useful 
structures such as microstrips on dielectric substrate used in monolithic microwave 
integrated circuits (MMIC’s) are not available (Ref. 2: pp. 256-262]. 

Knowledge of the transient signal behavior on microstrip transmission lines is es- 
sential for the design of MMIC’s at high switching speeds or high frequencies. High- 
speed time-domain measurements must be used to properly understand and model 
the transient response. The microwave techniques available for picosecond pulse prop- 
agation characterization are scattering parameter measurements, high-speed sam- 
pling oscilloscope measurements, and picosecond photoconductor measurements. The 
scattering or s-parameter measurement method is a small-signal frequency-domain 
technique which is widely used by microwave network analyzers. High-speed sam- 
pling oscilloscope measurements can characterize low-level signals and have a higher 
signal-to-noise (S/N) ratio than the s-parameter m»thod. Picosecond photoconduc- 
tor measurements are a new high-frequency measurement technique made possible by 
advances in laser technology. The generation of optical pulses with sub-picosecond 


duration are shorter than those that can be generated solely by electronic means [Ref. 











3: pp. 117-124]. An opto-electronic transducer such as a photoconductor permits the 
conversion of these sub-picosecond optical pulses to picosecond electrical signals. The 
significant advantage of this technique is the photoconductors can be integrated on 
the substrate material to facilitate very high-speed measurements with extremely high 
S/N and sensitivity to microvolt signal levels. The disadvantage of the s-parameter 
and sampling oscilloscope measurement methods is they suffer from poor connections 
in terms of high-frequency signal transmission on the substrate. 

There is an urgent need for the development of computer-aided design (CAD) 
models tailored specifically to the special demands of MMIC technology. MMIC’s are 
not readily tunable after manufacture, and the time and costs involved in a design 
are high. In addition to the development of new models for CAD design, a diagnostic 
testing procedure that permits on-wafer characterization of MMIC’s before dicing the 
wafer into individual chips is highly desirable [Ref. 4: p.14]. 

It will be useful to formulate a digital signal processing framework to solve this 
modeling problem. The motivation for this approach is simple. The availability of in- 
put and output time-domain measurements invites the transition from an equivalent 
circuit model to a parametric model. For our purposes, we will define a parametric 
model as a discrete mathematical description of the actual propagation mechanism of 
the system in terms of specific model parameters or coefficients. Three particularly 
important discrete-time parametric models are the moving-average (MA), the autore- 
gressive (AR), and their combination, the autoregressive moving-average (ARMA) 
model. These models are described as linear transformations of respective input- 
output time series. 

The model-building process is determined from a priori (structural) knowledge 
and a posteriori (measurement) knowledge of the system. The basic approach taken 
in this dissertation, with respect to the modeling of microwave integrated circuit 
(MIC) interconnections, is to represent them as linear time-invariant systems. This 
allows linear algorithms to be applied to the solution of the modeling and parameter 


estimation problem. The structura] knowledge of the microstrip line permits the use 


of empirical and analytical expressions to formulate lumped-distributed equivalent 








circuit models which can characterize the system. Rational network functions are 
evaluated from: the circuit models and discrete linear transformations are employed 
to produce the parametric models. 

Another main emphasis of this research is on the system identification aspect. 
Identification is defined by Zadeh (Ref. 5: pp. 856-865] as: “the determination on the 
basis of input and output, of a system within a specified class of systems, to which the 
system under test is equivalent.” Consequently, the identification problem is reduced 


to that of model parameter estimation based upon measurement knowledge. 


A. HISTORICAL BACKGROUND 

Numerous authors have focused their research on various computational 
approaches to analyze and model microstrip on different substrates [Refs. 6,7]. Ac- 
curate and complete analysis, however, requires elaborate mathematical models and 
time intensive numerical techniques. Quasi-TEM mode calculations combined with 
frequency-dependent expressions can achieve accuracies within one percent of the full- 
wave analysis [Ref. 8]. Computer-aided design (CAD) programs are available with 
the capability of both synthesis and analysis of microstrip lines. The extensive de- 
velopment of closed-form expressions for the frequency-dependent effective microstrip 
permittivity eeys(f), has produced rapid computation of the microstrip parameters. 
However, the empirical equations that have been derived for the frequency-dependent 
propagation terms do not lead to an equivalent circuit model with realizable series and 
shunt circuit impedances. It is these equivalent impedances that relate to the physi- 
cal microstrip properties. Therefore, the effective microstrip permittivity of different 
dielectric substrates is critical for the synthesis and validation of IC interconnections 
models. 

Recently, optical techniques have been successfully used in the characterization 
of microwave devices and integrated circuits. Frequency-domain measurements have 
been performed using electro-optic probing of a microstrip line. In this work, the 
microwave signal is launched onto the circuit using coplanar waveguide (CPW) con- 


tacting probes. Hung et al. [Ref. 9] have developed an on-wafer GaAs MMIC 











measurement system using a picosecond pulse sampling technique proposed by Aus- 


ton {Ref. 10]. This optoelectronic characterization has been demonstrated to achieve 
a broad-band frequency response for both the magnitude and the phase of a K,-band 
(27-40 GHz) MMIC. 

The problem of system identification is the determination of a mathematical 
model that provides an optimum characterization for a system or process. Usually 
our knowledge of the system is limited by observable input-output measurements. 
The inverse scattering problem as applied to electromagnetic systems has close con- 
nections to several signal processing concepts such as the design of digital filters, the 
development of linear prediction algorithms and their lattice filter implementations. 
In particular, transmission-line systems on planar dielectric materials can be inter- 
preted as layered, wave scattering structures. Two successful analogous applications 
have been in speech and seismic signal processing. In speech, the vocal tract has been 
modeled as an acoustic tube of varying cross-sectional area [Ref. 11]. The inverse, 
or model identification, problem is to determine the medium properties from its re- 
flection response measured at an observable boundary, to some incident input signal. 
This also relates to the seismic problem. Here, an impulsive input to the earth. such 
as an explosion, will generate seismic waves propagating downwards. Reflections are 
produced as the wave encounters respective earth layers. From the resulting reflected 
waves, the layered structure beneath the surface can be identified. 

Bruckstein and Kailath [Ref. 12] have specified conditions on the layered structure 
that permits the use of recursive parameter estimation algorithms. Their first algo- 
rithm takes scattering data and processes it to identify a unique layer and then, at each 
iteration, replaces the data by a set of “synthetic” scattering data. This procedure 
is called the layer-peeling method. The second approach, called the layer-adjoining 
method, propagates the original scattered data through the previously identified layer 
to determine information about future layers. The layer-adjoining method is the 
Levinson-Durbin alogrithm for solving the Toeplitz normal equations. For lossless 
structures, the layer-peeling method is the Schur algorithm. The primary difference 


between the layer-adjoining and the layer-peeling process is that the latter method 








has avoided the requirement to compute inner products. These algorithms will pro- 
vide background for the derivation of model parameter estimation algorithms that 
will be presented in Chapter V. 

Bruckstein, Levy and Kailath (Ref. 13] have discussed several classes of physical 
models that are equivalent descriptions of a lossless scattering media. They presented 
how a lumped circuit model for a uniform lossless transmission line, described by par- 
tial differential equations, can be formulated into a discretized wave scattering layered 
model. This model] illustrates how propagating waves through a lossless transmission 
line can be described as a multi-layered medium. Each layer is characterized by its 
impedance function which directly relates to a layer reflection coefficient. In Chap- 
ter IV, multiple microstrip discontinuities are compared to the multi-layered model: 


however, loss and dispersion are considered. 


B. OBJECTIVES OF THE RESEARCH 

The first objective of this research is to develop an equivalent circuit model for 
a dispersive lossy microstrip transmission line which is compatible with the standard 
circuit analysis and design techniques, including computer-aided design (CAD) tools. 
The proposed circuit model will then be extended to describe abrupt width dimen- 
sion changes, called impedance discontinuities, of the microstrip line. The transient 
analysis of these equivalent circuits will simulate the picosecond pulse propagation 
on dispersive lossy microstrip lines and the effects of impedance discontinuities can 
be modeled. The circuit models will be verified using photoconductor measurements 
from IC interconnections and experimental results from fabricated microstrip test 
structures. 

The second objective of the thesis is to derive ARMA and AR parametric mod- 
els for three microstrip test structures. These include 1) an impedance matched 
microstrip transmission line, 2) a cascaded microstrip step discontinuity, and 3) a 
multi-section microstrip step discontinuity. The impedance matched microstrip line 
describes the typical MIC interconnection [Ref. 3: pp. 18-19]. As discussed eariler, 


the formulation of these parametric models provides an opportunity to exploit existing 











parameter estimation algorithms, as well as develop new algorithms to characterize 
transient signal propagation on MIC’s. 

The last objective is the development of parameter estimation algorithms. This 
research work will focus on the estimation of both ARMA and AR model parameters 
from .he impulse response of an electrically short microstrip line section. Addition- 
ally, the system identification of the multi-section microstrip structure will also be 
investigated. The detection of each section impulse response will be complicated 
by pulse dispersion, propagation loss, and multiple reflections due to the impedance 
discontinuities. A new layer-probing algoritnm will be presented to overcome these 
difficulties. Computer simulations of the algorithms are performed to validate their 


performance. 


C. ORGANIZATION OF THE REPORT 

Chapter II derives a lumped-distributed equivalent circuit that includes the ef- 
fects of microstrip loss and dispersion. The bandwidth of the propagating pulse will 
establish a maximum wavelength of interest. The physical length of the IC intercon- 
nection being considered will be very much less than a quarter-wavelength. Under 
this condition the lumped-element circuit approximation is used to characterize the 
microstrip line section. 

Chapter III presents the parametric models for the impedance matched quarter- 
wavelength microstrip transmission line. An ARMA parametric model is derived 
from the network function of the equivalent circuit model. The ARMA model coef- 
ficients are shown to be directly related to the equivalent lumped capacitance and 
inductance of the circuit model. An AR parametric model will also be presented. 
Relationships are given to evaluate the effective microstrip permittivity directly from 
the ARMA/AR model parameters. 

In Chapter IV, microstrip discontinuities will be considered. Equivalent circuit 
models are derived for a cascaded microstrip step discontinuity. Similarly, ARMA and 
AR parametric models are developed using discrete transforinations of the equivalent 


circuit network functions. 











Parameter estimation algorithms are the emphasis in Chapter V. Several esti- 


mation algorithms are developed for both deterministic and stochastic data. Three 
deterministic based algorithms are presented. A weighted least squares (WLS) algo- 
rithm is derived to solve the ARMA model parameters from a finite-length impulse 
response. An alternative technique will analytically approximate a transfer function 
of the microstrip section using measured input-output rise and delay times of a tran- 
sient pulse. Finally, AR model parameters are estimated from the impulse response 
using the Schur algorithm [Ref. 14]. The stochastic algorithms will assume that a 
white noise source is applied to the input of the microstrip structure. The Schur 
algorithm will estimate the AR model parameters from the output data. When input 
and output random data are available, an ARMA parameter estimation algorithm 
based on a generalized Mullis-Roberts (M-R) criterion [Ref. 15] is employed. How- 
ever, if only the output data can be measured, a modified two-stage least squares 
algorithm is presented to estimate a second-order ARMA model. Finally, new sys- 
tem identification algorithms for the multi-section microstrip step discontinuity are 
presented. 

Chapter VI presents simulation and experimental results. Waveform compar- 
isons are made between measured data, equivalent circuit simulations, and paramet- 
ric model simulations for each microstrip test structure. The performance of the 
parameter estimation algorithms will be investigated. 

Chapter VII is a summary of the significant contributions presented in this dis- 
sertation. It draws conclusions from the results and proposes some important future 
directions for this research. 

Three appendices are included. Appendix A contains an alternate proof of the 
Schur algorithm. Appendix B contains the equivalent circuit model listings that 
are uscd by the transient analysis program. Appendix C contains listings of the 


FORTRAN programs used in the simulations presented in this report. 








II. MICROSTRIP CHARACTERIZATION 


In this chapter equivalent circuit models are derived for a lossy dispersive mi- 
crostrip transmission line. The propagation characteristics of the transmission line 
are modeled by lumped-distributed equivalent circuit models. In order to derive the 
expressions for the equivalent lumped circuit elements. the models are defined for 
a maximum frequency of interest which is restricted by the physical length of the 
microstrip line and the dielectric constant of the substrate. The empirical equations 
described in the literature for the propagation parameters do not usually lead to an 
equivalent circuit model [Ref. 16]. 

The equivalent circuit models presented are compatible with standard circuit 
analysis and design techniques including the use of computer-aided design tools such 
as PSPICE [Ref. 17]. The proposed models will include the effects of dispersion 


and loss at microwave frequencies. 


A. IC MICROSTRIP INTERCONNECTIONS 

Integrated circuit interconnections can be described by microstrip transmission 
lines because their geometries are similar. The abrupt dielectric interface shown 
by the open microstrip geometry in Figure 2.1 makes it incapable of supporting 
a single model of propagation. However, microstrip propagates the bulk of its 
energy in a field distribution which approximates the transverse electromagnetic 
(TEM) mode and is usually referred to as the quasi-TEM or quasi-static mode 
(Ref. 18]. Several computational approaches are available utilizing quasi-TEM mode 


calculations combined with closed-form frequency-dependent expressions {Ref. 6]. 


It is this latter approach which is considered in this work. 








The characteristic impedance of a TEM transmission line is described by 


ne ff an 


where L and C are the equivalent lumped inductance and capacitance of the mi- 
crostrip section, respectively. With microstrip geometries, the same type of di- 
electric substrate is used below the conductor. However, there is air above the 


conductor. 


Top Conducting Strip 






Dielectric Substrate 


Ground Plane (Conducting) 


Figure 2.1 Generalized open microstrip geometry 


The dielectric constant used in the design must take into consideration the dielectric 
constant of air (€- = 1) and that of the substrate material. The effective microstrip 
permittivity ¢e¢¢ will be .ntroduced. This quantity is unique to mixed-dielectric 
transmission line systems and it provides a useful link between different wavelengths, 


impedances, and propagating velocities. 





An expression for the static-TEM effective microstrip permittivity has been calcu- 


lated by Owens [Ref. 19]: 


—0.555 
é¢>¢+1 6-1 h 
tgp eS 2.2 
where ¢, is the relative dielectric constant of the substrate. Using the effective mi- 
crostrip permittivity, the characteristic impedance at TEM frequencies is calculated 


by 





60 Sh w 
— se a= ‘ ): 
Zo _ In (= + x) u<h (2.30) 
1207 u 





a (= + 1.393 + 0.667 In (- rs 1444) roe (2.30) 
where w and h are microstrip conductor width and substrate thickness, respectively 
(Ref. 20]. A correction factor is applied to account for the fringing fields associated 
with a finite conductor thickness. Bahl et al. (Ref. 21] have introduced an effective 
width parameter into (2.3) in order to improve the Zp calculation. The effective 


width we, replaces the w in equation (2.3) 


1.25t 4nrw Uw ] 
We = w+ —{14I1n : = — (2.4a) 
T t h Qn 


1.251 2! ’ 
we = w+—(14+In zal : = = (2.45) 
T t h Qn 


where ¢ is the conductor thickness. 
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B. FREQUENCY-DEPENDENT MICROSTRIP MODELS 

The quasi-TEM analysis of microstrip transmission lines mentioned above he- 
gins to lose accuracy at high microwave frequencies. The characteristic impedance 
and the effective microstrip permittivity are dependent upon the dielectric thickness- 
to-guide wavelength ratio, 2th/A, [Ref. 22]. In addition, different loss mechanisms 
become important and the attenuation function is also frequency dependent. This 


frequency-dependency of the microstrip is caused by a hybrid mode of propagation 
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that describes a coupled version of both transverse electric (TE) and transverse 
magnetic (TM) modes. 
The propagation group velocity of a signal also depends on the frequency- 


dependent effective microstrip permittivity as 


een as 9.5 
Ug(f) eG) (2.5) 


where c is the velocity of light in a vacuum. Pulses which have a spectral component 
above the TEM mode frequency regime will be dispersed because the higher har- 
monics of the pulse will travel at a slower phase velocity than the lower harmonics 
[Ref. 18]. Therefore, the phase constant (;3) is a non-linear function of frequency. 
which leads to the phenomena of dispersion. 

Several methods for evaluating the frequency-dependent effective microstrip per- 
muttivity are available [Refs. 23. 24]. However. Yamashita et al. [Ref. 7] have 
derived an expression by curve fitting the datu obtained from a full-wave analysis 


as 
eps (f) = Jeege t (Va - veg) / (1 +4F7?) (2.6) 
where 


h 


If the substrate material is lossless, the phase constant becomes 


p= (£)nve=n [54 (142 toe (14+ Z)) 


omy SeeTP 7 
js (2.7) 
c 
However, finite resistivity of the conductor and finite conductivity of the substrate 
introduce attenuation. At low microwave frequencies the conductor loss factor a, 
is given by [Ref.18: p. 90] 


[ 0.072/f Ag 
= w Zo 


ac 


dB/A, (2.8) 
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where f is in gigahertz. For higher frequencies, Puce. et al. [Ref. 25] have calculated 


the conductor loss as 











c= mil CG) Ita tye tah eae 
for " < - 
Ac = oS [1 = (=) + a + (in a _ ~)]: (2.96) 
for = < - 2 
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where a, is in dB/cm and 





w =w+Aw 
t 4rw 2t w il 
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and R, is the surface skin resistance in ohms given by 
eee oni (2.94) 
oc 


where yz: is the conductor permeability and o, is the conductor conductivity. An 


expression for dielectric loss ag has been derived by Hammerstad and Bekkadal 

[Ref. 26] as 

er(€ege¢ — 1) tan d 
ef flO = 1) 





ag = 27.3 dB/q (2.10) 





where tan 6 is the lcss tangent of the substrate. For lossy semiconductor substrates. 
the general attenuation and phase constants for the microstrip are derived froin an 
analysis of a transmission line section modeled by series and shunt resistances to 


represent the conductor and dielectric loss mechanisms [Ref. 27] as 


(32+ 4a?) (63 + 4a’) — Bs + 4a,.aq 

2 3 

(32 + 4a2) (a ze 4a3) + Be — 4a,aq 
ee 


C. PICOSECOND PULSE PROPAGATION 

One of the research objectives is to accurately characterize picosecond pulses 
propagating on microwaye IC interconnections. Researchers have studied the ef- 
fects of dispersion on picosecond pulses by performing computer simulations in the 
frequency domain [Refs. 28,29]. An input time-domain signal v(t, 0) is transformed 


using Fast Fourier transform (FFT) techniques to produce the input frequency spec- 


trum, F{ v(t.0) \. A frequency-dependent propagation function +(f) acts upon each 


spectral component as 


F{u(t.0) be WW) l 


where the propagation function is 4(f) = a(f)+ j3(f) and / is the propagation 
distance along the microstrip. The distorted temporal pulse is determined using an 


inverse Fourier transform as 
v(t.) = Fl F{o(t.0) bela) + j3(f)] \. 


The approach taken in this work is similar, but the exponential propagation term 
is replaced by a discrete parametric model of the microstrip section. Specifically, 
the phase term 6(f) is modeled by a digital transfer function H(z), while a discrete 
loss multiplier describes the attenuation constant a(f) of the microstrip line section 
at a maximum frequency of interest. Chapter III will introduce the parametric 


models for several microstrip transmission line structures. 





D. DERIVATION OF EQUIVALENT CIRCUIT MODELS 
A microstrip transmission line is approximated by cascading several electrically 
short sections [Ref. 30]. A single section of microstrip line having a characteristic 
impedance Zg and terminated in a matched load impedance is shown in Figure 
2.2a. The load impedance Zn4}), is the equivalent input impedance of the (n + 1)- 
th section. The normalized input impedance to the n-th section is derived (Ref. 30] 
as 
Zn - Zn+l cosh(y!) + Zo sinh(+/) (2.12) 
Zo = Zn41 sinh(yl) + Zo cosh(y/) 
where the propagation function, y = a + j8, describes the attenuation a, and phase 


shift @ of the n-th section. Equation (2.12) is rewritten as 


2 Aas + tanh(7/) 
Zz = — (2.13) 





Zo tanh(7/) 


The hyperbolic function, tanh(7/), is expanded to 


tanh(al) + j tan(l) 








=. 2.14 
ean) 1 + j tanh(a@l) tan(8!) ete) 
Substituting (2.14) into (2.13) yields 
f 
Z (4+) + tanh(al) 
$= a (2.15a) 
: 1+ ( x21 ) tanh(al) 
Zo 
where Zz 
n+l] : 
+ jtan(@l) 
Zn+1 Zo 
= 2.155 
(4+) (2.188) 


l+j (4+) tan(l) 
Zo 
A physical interpretation of (2.15) is illustrated in Figure 2.2b. The original n- 


th section is equivalent to three sub-sections of the same characteristic impedance 


cascaded together. The first and third sub-sections provides attenuation but no 
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phase shift, and are modeled by purely resistive attenuator pads. The total loss of 
the microstrip line is divided equally between each z-attenuator. The center sub- 
section has a pure imaginary propagation function and no attenuation. This type 
of reactive network corresponds to a lossless transmission line having a phase shift 
equal to the original section. 

In classical filter theory the reactive lossless, two-port network has been de- 
scribed as a constant-k lowpass filter. Two variations of this filter network will now 


be presented. A reactive m-network is used to derive the lumped-parameter 











(b) 


Figure 2.2 Microstrip line section equivalent circuit 


equivalent circuit model and a reactive T-network is used to derive a lumped- 


distributed equivalent circuit model. 
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1. Lumped-Parameter Equivalent Reactive I]-Network 
The lossless reactive network of Figure 2.2b is modeled by an equivalent 


m-network. The complex impedances are 


Z\ = ee (2.16a) 
) tan (5) 
J 2 

Z2 = jZo sin( Bl). (2.16) 


Equations (2.16a) and (2.16b) represent a capacitive and an inductive lumped re- 
active impedance, respectively. Solving for the susceptive (Bc) and reactive (Xz) 


lumped elements yields 


1 Bl 
Bo= > — 217 
C gven ( 5 ) (2.17a) 
Xp = Zosin( sl). (2.176) 


Assuming a very short section leugth (1 << 4,/4), the equivalent shunt capacitance 


and series inductance can be approximated as follows [Ref. 18] 


l 


a 2.18 
Cx 20 f Ay (2.18a) 
ea 
Le SS . 2.186 
ae (2.188) 


The guide wavelength , is expressed in terms of the frequency-dependent microstrip 
permittivity eess(f) as 
(2.19) 


C 
a 

pi V Ceff (f ) 

where f is the highest frequency of interest of the propagating signal [Ref. 18: 
p. 70]. The resultant lumped-parameter equivalent circuit for the microstrip line 


section is shown in Figure 2.3. 
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2. Lumped-Distributed Equivalent Reactive T-Network 
The reactive 7-network of Figure 2.3 can be replaced by an equivalent T- 
network as shown in Figure 2.4a. The series inductive and the shunt capacitive 


impedances are then given by 


l 
Z, = jZo tan (5) (2.20a) 
Zo 
= 2.206 
| Fsin( 3) mor 
Solving for the reactive and susceptive lumped elements of Figure 2.4 yields 
l 
X1, = 2Zotan ($=) (2.214) 
Xg 
1 2x 1 
= sin | — }. 22 
Be Z sn ( Vy ) (2.216) 





Figure 2.3 Lumped-parameter equivalent circuit model 


Assuming a section length that is less than a quarter wavelength, the lumped ele- 


ments can be approximated as follows [Ref. 18]: 
, Zol 


Lrs — 2.224 
rey i ( ) 

l 
Cr xs ——-_. 2.226 
T ZoPhe ( ) 


Next an equivalency between the series lumped inductance (Ly) and a lossless 


transmissiom line section is derived. Referring to Figure 2.4b. we define equivalent 


17 





impedances, Zeg, with respect to the shunt capacitance C’r. The load impedance 
Zr is equivalent to the characteristic impedance Zp. Therefore, the equivalent 
impedance yields 


Zeq = ZL + jwlLr. (2.23a) 


Using (2.15), the n-th section input impedance for a lossless (a = 0), high impedance 


transmission-line section is 


Zn+1 + jZo tan(Pl) 





A= FZ (2.23b) 
1+ j— tan(8l) 
Zo 
Assuming |91| < 1, the input impedance can be approximated by 
Zn & Zn+1 + jZo(8l). (2.23c) 


Equation (2.23c) is now equated to the equivalent input impedance Ze,. 





Figure 2.4(a) Equivalent T-network equivalent circuit 


Assuming a small section length and taking the imaginary terms in (2.23a) and 


(2.23c) into account yields 


Lr Zo (=) =Zota (2.24) 


where tg is the propagation delay through the lossless transmission-line section. 
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Each lumped inductance Ly can now be replaced by an equivalent distributed trans- 


mission line as shown in Figure 2.5. 





Figure 2.5 Lumped-distributed equivalent circuit model 
The accuracy of the circuit model is dependent on the microstrip line length 
I, and the maximum frequency bandwidth wmaz of the propagating transient sig- 
nal. The equivalent lumped-distributed circuit presented in this section is valid for 
electrically short, high impedance transmission lines which satisfy the conditions 
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E. ESTIMATION OF CHARACTERISTIC IMPEDANCE 
Picosecond pulses propagating on the IC interconnections are described as finite 


energy signals and their maximum cumulative energy is defined by [Ref. 31: p. 35] 


+00 
Emaz,= >. \v*[l,n][’ (2.25) 


n=-c0 


where vt [I,n] is the sampled right-propagating voltage measured at the I-th length 
along the microstrip. To discuss pulse propagation on a lossy microstrip line, we first 
obtain the transmission-line equations by applying Kirchoff’s voltage and current 
laws to the alternative distributed equivalent circuit as shown in Figure 2.6. When 
low-loss, high frequency conditions are assumed, the general solution for the voltage 
and current at the /-th length is solved as [Ref. 32: pp. 437-445] 


v(1,t) = vt (: - =) +07 (: + =) (2.262) 
i(l,t) = oa jo (: - ~) +7 (: + ~)| (2.260) 


where the superscripts (+) and (-) denote the right- and left-propagating voltage 
waves, respectively. When we compare these equations, the right-propagating cur- 


rent and voltage are given by 


1 i 
+ Ss a te 
2 (i,t) = Za” (: =) (2.27a) 
vt (1, t) = vt (: ~ ~) (2.276) 
U% 


where vg is the group velocity. The instantaneous power associated with the right- 


propagating wave is 


p* (i,t) = et(1,t) i*(Lt) (2.28a) 
featecak\\ 
24 all (2.288) 
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Since the pulse rise time is very fast (< 10 ps), the derivative of the cumulative 
energy curve (i.e. the instantaneous power) is approximately constant. Therefore, 


the maximum power of a causal, transient pulse is given by 
1 a 2 
Ping os Wo let (1, n]j (2.29) 


where N is the number of voltage sample values. By equating (2.28b) and (2.29), 


we have the characteristic impedance given by 


2 
= jut [l, fil peak | 


Zo = 2.30 
ay eee) 


where v*[I,n]peak is the peak amplitude voltage sample of the propagating pulse. 


i(f,1) RCAL) LAL) i(f£+Al,t) 







v(f,t) v(f+Al,t) 


a 


( Lumped-elements expressed as per unit length ) 


Figure 2.6 Lossy transmission line equivalent circuit (After Ref. 32) 


In summary, distributed-lumped equivalent circuit models were derived for an 
electrically short length of lossy microstrip transmission line. The lumped circuit 
elements are evaluated at the maximum frequency of a bandlimited input signal. 
The validity of the equivalent circuit models are controlled by the physical length 
of the microstrip line. Therefore, we are interested in picosecond pulse widths 
and specific IC interconnections lengths which satisfy the lumped-element modeling 
criterion. Chapter III will use the proposed circuit models as a foundation to derive 


the parametric models for an IC interconnection. 


Ill. PARAMETRIC MODELING 


The rationale for using a parametric modeling approach to approximate a mi- 
crostrip transmission line is intimately related to the identification of the physical 
microstrip properties. This chapter will present both autoregressive moving-average 
(ARMA) and autoregressive (AR) models of the equivalent reactive networks pre- 
viously derived for the IC interconnection. A second objective is to show how 
the model parameters are related to the effective microstrip permittivity e-sy, the 
characteristic impedance Z, and the propagation group delay ry cf the microstrip 


structul. 


A. RATIONAL TRANSFER FUNCTION MODELS 
1. ARMA Parametric Model 
In this model, the discrete-time input sequence z{n], and the microstrip sec- 


tion output sequence y[n], are related by the linear difference equation [Ref. 31} 


P 


q 
y{r] = So bea[n — k] — So axyln — fi. (3.1) 
k=0 


k=1 


where 5, are the moving-average (MA) model parameters and a, are the autore- 
gressive (AR) model parameters. The flowgraph representation of (3.1) is shown in 
Figure 3.1. This difference equation describes an ARMA model of order (p,q) and 


the corresponding system transfer function is 





Se (3.2) 
1— Yo agz-* 
k=] 
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2. AR Parametric Model 


When all the 0, coefficients, except bp) = 1, are zero in the ARMA model, 
equation (3.1) reduces to 


Pp 
y[n] = z[n] — D> axy[n — &] (3.3) 
k=1 


and the corresponding AR or all-pole transfer function is given by (Ref. 33: p.111] 


1 1 
1— Sayz-* . 
k=] 
INPUT SEQUENCE 


x{n] © 






OUTPUT SEQUENCE 
y(n] 


Figure 3.1 Realization of the ARMA/AR parametric models 


In Figure 3.1, a flowgraph realization of this model is shown by the weighted feed- 
back signal path. 











B. IC INTERCONNECTION PARAMETRIC MODELS 
1. Derivation of the ARMA Model 
The network transmission function for the reactive 7-network of Figure 2.3 


yields 


w9 


H(s) = —=————_ 
(s) s? + 2wos + w? 


(3.5) 
where wo = (Lae) and ¢ = (2ZoCxw)7!. This network function describes a 
bandlimited frequency response for a second-order lowpass filter. Here we use the 
impulse invariant design method to transform (3.5) into a digital transfer function 


(Ref. 34]. Equation (3.5) can be expressed in the partial-fraction expansion form as 





2 Cy 
H(s)= 0 —5 (3.6) 
k=] k 
where 
Ch «0 ett /2 


and the complex poles are 


dy,2 = —Cwo + jun V1 — ¢?. 


Taking the inverse Laplace transform, the corresponding unit impulse response be- 


comes 


2 
h(t) = S>Cyet#t for 120 (3.7a) 
k=1 


and its discrete representation is 


2 
h[n} = D> Cedar (3.76) 
k=1 


where T is the sampling interval. 


24 











It should be observed that for high sampling rates (say, T = 1.25 ps) the digital 


filter has an extremely high gain. For this reason (3.7b) is expressed as 


2 
h[n] = TS Cpe? (3.8) 
k=] 


The digital transfer function is then obtained by taking the z-transform of (3.8) as: 


2 


7 TC; 
Harma(z) = do PST VET (3.9) 
After rearranging the terms and simplifying, we have 
by27} 
Harmalz) (3.10) 


~ 1 ~ayz7! — agz7? 


where the filter coefficients, also called the ARMA model parameters, are: 


_ wl  _fuy T 7 
b} = Were (wo cos G —up V1 — ¢? i) (3.11a) 
2 e—Swol cog (“ovi —¢? T) (3.118) 


ay = —e2suo 7. (3.11c) 


a} 


The filter coefficient 6; is directly linked to the AR filter coefficients a; and ap. 
Given a2 and the sampling period T, (3.11c) yields (uy. Substituting this result 
into (3.11b), a value for wo fi-@ can be found. Finally, specific values for wo and 
¢ are easily determined and the numerator filter coefficient 5; can be calculated. 
Therefore, the MA model parameter in (3.10) can be obtained from the AR model 
parameter estimates. In Chapter IV we will take advantage of this computational 


result in our discussion of parameter estimation algorithms. 
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2. Derivation of the AR Model 


The AR parametric model is derived from the lumped-distributed equivalent 


circuit model of Figure 2.5. A network reflection function I'(s) is solved as [Ref. 35] 


r(s) = Yo(s) — ¥z(s) 


= Ya(s) 4 ¥i(s) a 


where the equivalent admittances are given by 


Yo(s) = z 
Yz(s) = Yc(s) + Yo(s) 
1 
=sCr+ 


Z 
After simpilfying and rearranging (3.12), we have 
(2) 
ae or : (3.13) 
l+s (S) Zo 
The network transmission function is then defined as 


T(s) =14+T(s). (3.14) 


By substituting (3.13) into (3.14), a first-order, single time constant transfer func- 
tion is achieved: 


T(s) = (3.15) 





v4 
where 6 = (ae G ). Taking the inverse Laplace transform yields the unit impulse 
oO™ 


response 
har(t) =e 6! for t>0. (3.16) 


The z-transform of the corresponding normalized sampled transmission impulse 


response is given by 


1 
Har(2) = ——z—— for [z| > le~7>|. (3.17) 
l-—e “"z 
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where T,, is the normalized sampling interval. In order to satisfy the Nyquist sam- 
pling criterion, the sampling intervals of the normalized and denormalized trans- 


mission responses must be related by [Ref. 36: pp. 341-342] 


7 = eT (3.18) 


To adjust the dc gain of (3.17) we can multiply the digital transfer function H4p(z) 
by the ratio of (3.15) evaluated at s = 0 to (3.17) evaluated at z = 1 ( which is the 
dc value). Thus, the AR transfer function, after being adjusted for the dc gain, 
becomes 


] -Tn 


as for |z| > |e~Z*|. (3.19) 


Har(z) = ——= > 
J—e tn 2-1 


3. Effective Microstrip Permittivity 
The effective microstrip permittivity is directly related to the ARMA model 
parameter a2, the sampling interval T, and the physical microstrip length /. Using 
(3.11c), the propagation group delay is solved as [Ref. 18: p. 3] 





(3.20a) 


= Tey (3.206) 


where the group delay of the reactive network is Ty = os ae aes The propagation 
group delay can also be related to the effective microstrip permittivity by [Ref. 18: 
p. 62] 

Veeif (3.21) 


%= é 
. c 





Substituting (3.21) into (3.20b) yields the effective microstrip permittivity 


ae |. (3.22) 


“eff = | naa) (0) 











We now derive an expression for the effective microstrip permittivity from the 
AR model parameter. From (3.19), we recognize the first-order AR filter coefficient 
iS 


aze ln, (3.23) 


The effective microstrip permittivity is now solved in terms of this AR model pa- 


rameter and the propagation group delay: 


2T 
= Te (3.246) 


Substituting (3.21) into (3.24b) yields the effective microstrip permittivity 


2 
€eff = mail ; (3,25) 
4. Modeling Propagation Loss 

Both the ARMA and AR digital transfer functions were derived from a loss- 
less equivalent circuit. In order to properly characterize a picosecond pulse prop- 
agating on the IC interconnection, we must include a loss mechanism in the para- 
metric models. The total propagation loss (in dB) experienced by a high-speed 

transient pulse will be defined by [Ref. 37: p. 7] 


Lap = 10logy, lore | (3.26) 
where a; is the total attenuation factor (in nepers/m) and /| is the propagation 


distance (in meters) along the microstrip. The attenuation factor can be solved by 


combining the theoretical conductor and dielectric attenuation factors from (2.8b) 


and (2.10), respectively. 















In summary, this chapter has presented ARMA and AR parametric models for a 
lossy dispersive microstrip transmission line. Analytical expressions were derived to 
relate the effective microstrip permittivity to the model parameters at a maximum 
frequency of interest. We have assumed no discontinuities on the IC interconnection 
because of impedance matching. In Chapter IV, we will extend the equivalent circuit 
and parametric models for a cascaded microstrip step discontinuity. Chapter V will 
present several parameter estimation algorithms which will be used to approximate 


the effective microstrip permittivity of the IC interconnection. 








IV. MICROSTRIP DISCONTINUITY MODELING 


Numerous approaches have been made to use equivalent circuits tc model the 
microstrip discontinuities. Whinnery and Jamieson [Ref. 38] used Hahn’s method 
to match electromagnetic-wave solutions across discontinuities, Oliner [Ref. 39] 
used Babinet’s principle to describe stripline discontinuities, and Menzel and Wolff 
[Ref. 40} calculated the frequency-dependent properties of various microstrip dis- 
continuities. 

A discrete parametric method is presented in this chapter to model the reflec- 
tion and transmission characteristics of a cascaded microstrip step discontinuity, 
assuming the quasi-TEM mode of propagation. Equivalent circuits are developed 
using the shunt capacitance circuit model presented by Gupta and Gopinath [Ref. 
41). Both reflection and transmission network functions are derived from the circuit 
models. ARMA and AR digital filters are described using discrete transformations 
of the respective network functions. This method has the advantage that complex 
microstrip circuits containing discontinuities can now be modeled by linear differ- 


ence equations. 


A. MICROSTRIP DISCONTINUITY EQUIVALENT CIRCUITS 
1. Single Step Discontinuity 
A single microstrip step discontinuity exists at the junction of two microstrip 


lines having different widths and characteristic impedances. This type of discon- 


tinuity is seen in the design of microwave matching transformers, couplers, filters 


and transitions (Ref. 41]. 











The geometric configuration of the step discontinuity and its equivalent cir- 


cuit are shown in Figure 4.1. 


—~——O 


Figure 4.1 Microstrip step discontinuity and the equivalent circuit 
The effect of the total discontinuity inductance (Lg) may be separated into L, and 


Iz as 
Lw1 
Ll) = ——_——_L H 4,] 
a ca ae (H) (4.1) 
Lw2 
Lo = ————-L H 4.1b 
ak ea ER (H) (4.18) 


where Lyi and Ly2 are the inductances per unit length for the microstrip of width 


W, and W?2 given by 


CO raya 
Lum = Vin, m=1,2 (H/m) (4.2) 





where Zom is the characteristic impedance and €,f,,, is effective permittivity of the 
width Wy, and c = 3 x 108 m/s. The closed-form expressions for Cg and Lg have 


been derived from curve fitting numerical results. These expressions are [Ref. 41] 


Ca 





= (10.1 loge, + 2.33) — 12.6 loge, —3.17 (pF/m) 
2 


: 


(for €- < 10; 1.5 < Wi/W2 < 3.5) (4.3a) 
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La _ Ww, WwW, W, 
Tee 40.5 GE —1. 0) - 75 log 7 + 0.2 Ca —1. 0) (nH/m) (4.36) 


where e€, is the relative dielectric constant and h is the substrate thickness. Equation 
(4.3a) yields a percentage of error less than 10 percent and equation (4.3b) has an 
error less than 5 percent for W)/W2 <5 and W2/h = 1.0. 


2. Cascaded Step Discontinuity 

The cascaded microstrip step discontinuity is formed by combining two sin- 
gle step discontinuities as shown in Figure 4.2(a). Short (< 4,/4) lengths of high 
impedance (narrow width) microstrip line will behave predominantly as a series in- 
ductance (L,). Whereas, a very short < A,/4 length of low impedance (wide width) 
line will act predominantly as a shunt capacitance (C,) [Ref. 18: pp. 212-216]. 
End-inductances (Z.) are also introduced in the low impedance line as its length 
approaches a quarter-wavelength. The predominant series inductive reactance of 

the line of length /2 is 
wl, = Zo2 sinh G <2). (4.4) 

Ag2 


Assuming the small-angle approximation such that (2712/2 < 7/4), equation (4.4) 


is approximated as 


wm 202 ly 


Ls 
aie 2 Ag2 





(4.5) 


where the frequency-dependent microstrip wavelength is 


c 


AQ5= Oooo—=—==. 
ae ees s(f) 


When high impedance widths (W2) are equal, the equivalent circuit can be described 
by Figure 4.2(b). 











The predominant shunt capacitive suspectance of the line of length 1) is 


Pee (= ) (4.6) 
: 1 











(a) 
L L 
L pa L L 
L | Ly 1 2 2 ; 2 L. 


STEP DISCONTINUITY STEP DISCONTINUITY 
(b) 
Figure 4.2 Cascaded Step Discontinuity and the Equivalent Circuit 


and for 271; /X91 < 7/4, the shunt capacitance is approximated by 


fT 
C,2-— 4.7 
f Zo1 Agi a0) 


Similarly, the end-inductances are obtained as 





Zoi hh 
Le& : 4.8 








The relatively low characteristic impedance Zo will cause a very small inductive 
effect (Ly1) at each step discontinuity. Therefore, at frequencies up to a few giga- 
hertz, the inductance L; in the equivalent circuit can be neglected [Ref.18: p. 217]. 
A modified equivalent circuit is formed by converting the low impedance equivalent 
T-network into a m-network as shown in Figure 4.3(a). The shunt capacitances 
C,/2 and the discontinuity capacitances Cy are added to form C,. Also, each step 
discontinuity inductance L2 is combined with their adjacent series inductance L, to 
form a new inductance L,. Assuming very short microstrip lengths and from (2.24), 
the new inductance can be approximated as a lossless, distributed transmission line 


section having a characteristic impedance of Zo2 and a time delay approximated by 





(4.9) 


Figures 4.3(b) and 4.3(c) show the complete lumped-distributed equivalent circuits 


for the cascaded microstrip step discontinuity. 


(b) 


Figure 4.3 Modified cascaded step discontinuity equivalent circuit 





(c) 


Figure 4.3 Modified cascaded step discontinuity equivalent circuit (continued) 


3. Capacitive Shunt Equivalent Circuit Model 
An alternative equivalent circuit is presented that is designed to yield an “all- 
pole” transmission network function. The equivalent z-network of Figure 4.3(b) is 
converted into a T-network and the series inductances are combined as shown in 
Figure 4.4(a). The combined inductances are approximated as lossless, distributed 


transmission line sections having time delays of 
fe (2L, + Le) 


t oon (4.10) 


The shunt capacitance equivalent circuit model is shown in Figure 4.4(b) where t) 


represents the propagation delay time. 


L. L. 
L —_ nt 


(a) 


Figure 4.4 Capacitive shunt equivalent circuit 
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(b) 
Figure 4.4 Capacitive shunt equivalent circuit (continued) 


4. Modeling Propagation Loss 
Chapter II presented a resistive t-network attenuator to model the total 
propagation loss in a single microstrip line section. This scheme will again be used 
for the equivalent circuits of Figures 4.3(c) and 4.4(b). Figure 4.5 illustrates how the 
half-symmetrical losses of the cascaded microstrip step discontinuity are modeled. 


Using (2.8b) and (2.10), a total propagation loss can be obtained for the micros:rip. 


Lossless Reactive 
Equivalent Network 





Figure 4.5 Lossy cascaded microstrip step equivalent circuit 


Next, we equally represent the total loss by two identical resistive 7-networks as 


shown in Figure 4.5. 
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B. PARAMETRIC MODELS FOR MICROSTRIP DISCONTINUITIES 
1. Reflection and Transmission Network Functions 
The characteristic impedance mismatch caused by the abrupt change in the 
width dimension will produce both reflected and transmitted traveling voltage wave- 
forms. The objective of this section is to derive reflection and transmission network 
functions which accurately characterize the cascaded microstrip step discontinuity 
propagation effects. Initially, a third-order transmission network function is derived 
from the equivalent circuit shown in Figure 4.3(c). This network function is trans- 
formed into an ARMA model. Next, a first-order approximation is solved from the 
single capacitive shunt equivalent circuit of Figure 4.4(b). This single-impedance, 
single-capacitor network function is then transformed into an AR model. 
a. Third-Order Transmission Function . 
The reflection network function for the 7-network of Figure 4.3(c), ex- 


pressed in the complex s-domain, is defined as [Ref. 38: p.217] 


_ Yoo(s) — ¥z(s) 
I'(s) = Yoo(s) + Yz(s) (4.11) 


where the admittances Yo2(s) and Yz,(s) are obtained as indicated in Figure 4.6. 


02 





Figure 4.6 Characteristic and load admittances 
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Substituting for ¥o2(s) and ¥z,(s) in (4.11) and manipulating the resultant expression 
yields 


(C2L.Z2y)s? + (2CyZey + Le)s 


T =i Sa eo 


(4.12) 


where the subscript (3) denotes a third-order network function. The transmission 


function is defined as [Ref. 38: p. 217] 


T3(s) =1+ T'3(s). (4.13) 


Substituting (4.12) into (4.13) yields 


(azn) * + (aap) #* (cree) 
2 ICxZ2,+ 5 2 
3 2 x“02 e 
ost (=) oe ( C2L.Z3, )s as (azz) 


b. First-Order Reflection and Transmission Functions 


T3(s) = (4.14) 





A first-order approximation of the reflection network function is deter- 
mined by solving (4.11) for the equivalent circuit of Figure 4.4(b). Evaluating 
(4.11) yields 


s (C,Zo2) 
s) = -—-————_—_——_ l 
NGS <Te 5 (CeZon) (4 5a) 
and the transmission network function becomes 
1 
T;(s) = (4.155) 


l+s (CxZo2) 


2. ARMA/AR Models for the Cascaded Step Discontinuity 
a. Autoregressive Moving-Average Parametric Model 
The first step is to perform a partial fraction expansion of equation (4.14). 
The s-domain expansion is generalized as 


ky K2e?® Kye~79 


ne Ged) Gam Goa 


(4.16) 
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where the three poles are: 


1 
d= -=——- 
: CxZo2 
d2 =-d + ju 
d, = —o0 — jw. 


Using the impulse-invariant design method, a digital filter transfer function is ob- 
tained from (4.16) as 
ky 
i42) = arr ee + 


2K2 cos(@) — 2K2e~7 T 2-1 cos(@ —w T) 


4.17 
1 — 2e-@ T cos(w T)z~! + e~ 41 Tz-?, ont) 


After considerable complex algebra the generalized ARMA transfer function be- 


comes a4 A 3 
Ts(z) = bo + bz" + bez 


= 4.18 
1 — ayz7! — agz7? — a3z ( ) 


=a 
where 

by = Ky + 2K c0s(8) 

ieee [Roe cos(6) + e7T (Ky cos(wT) + K2c0s(8 — wT))] 


by = eT [Ks +2Kye~77 cos(@ — wT] 


2e~77T cos(wT) + e~ UT 


ay= 
a3 = eT, 


The most significant filter coefficient a3 identifies the time constant (Zo2C;) of the 


network. In terms of the sampling interval T, this filter coefficient is 


2r 
a3 = exp (-zz5): (4.19) 








A reflection digital transfer function is obtained from the bilinear transformation of 


equation (4.15a) as [Ref. 34: pp.206-211] 





Ti(z) =Ti(s) 9/1) —27} 
= T (; + =| 
2C,Z2(1 - 27) 
(T + 2C,Zo2) + (T — 2CxZ2)=-} 


Ii(z) = “oe =) =) (4.20) 


T(z) =— 


and the transmission transfer function is evaluated as 


5 T(1+ 27) 
T(z) a (T + 2C,Zo02) + (T _ 2C, Z02) zal 
bo (1 + 27? 


which is a lst-order ARMA model. The parameters are defined as follows: 
a ) 
0 AP + 20,202 J’ 
yp = (_2CxZo2 
0 AT + 20,202 J’ 


pee T — 2C,Z02 
PN TLIC Tian 


b. Autoregressive Parametric Model 
A first-order AR digital filter of the transmission network function is ob- 
tained using the matched z-transform design method. The single pole of equation 
(4.15b) is directly related to the pole of the AR transfer function [Ref. 42: pp. 


666-667] as 
( 7 ) 
—_ ACrZo2) (4.21) 


~ (az) 
t—e \CrZo2/,-} 
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Tar(z) = 








c. Realization of the Lossy Parametric Model 
The proposed ARMA/AR reflection and transmission transfer functions 
were derived from lossless reactive networks. The propagation losses which are as- 
sociated with both the low and high impedance microstrip lines must be included in 
the parametric models. These losses are represented as constant “loss” multipliers. 
The value of the multiplier is defined by 


L = exp [—-(a J)] (4.22) 


where | represents the microstrip line length and a; is the total attenuation factor 
in nepers per unit length. Figure 4.7 describes the lossy parametric model for the 


cascaded microstrip step discontinuity. 


3. Parametric Models for Multiple Discontinuities 

In this section a parametric model is developed for a multi-section microstrip 
structure described by a combination of cascaded step discontinuities. Wave prop- 
agation in the multi-section microstrip structure exhibits several of the physical 
characteristics of the multi-layered earth model developed in geophysical signal 
processing [Ref. 35: pp. 306-315]. Similarly, two different propagating waves exist 
within each microstrip section, an incident (right-propagating) wave and a reflected 
(left-propagating) wave. Unlike the lossless multi-layered earth model, propagating 
signals in microstrip are influenced by frequency-dependent losses and dispersion. 
Energy conservation is established across each discontinuity boundary. However, 
energy loss will occur between adjacent discontinuities. A linear time-invariant 
system is used to describe each discontinuity boundary. At each discontinuity an 
incident signal is partially reflected and partially transmitted into the next section. 
A proportional response of each component is given by the concurrent convolution 
of the incident signal with both a reflected and a transmitted impulse response of 


the linear system. A lattice-like system flowgraph relating right-propagating v7,(z), 
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and left-propagating v;,(z) signals at the m-th discontinuity boundary is shown in 


Figure 4.8. The reflection transfer function I',(z) 








REFLECTED RESPONSE 


rf. 


x; [ n ] 
INCIDENT SIGNAL 


Tfn] 


TRANSMITTED RESPONSE 


Figure 4.7 1st-order ARMA model realization of the cascaded step discontinuity 


and the transmission transfer function T,,(z) are the 1st-order ARMA models of 
(4.20a) and (4.20b), respectively. The overall parametric model for multiple dis- 
continuities may then be constructed by using the lattice as building blocks, and 
interconnecting the lattice sections with lossy delay lines. Figure 4.9 describes two 


adjacent discontinuity sections of a multi-section structure. Here an incident voltage 
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sequence oft Yin] initially enters a lossy delay line Dm(z), and a reflected voltage 


sequence oly Yin] is the resultant outcome after multiple reflections from the other 


discontinuities. 
Vv (+) (z) Vv (+) (z) 
m m +1 
(-) vi 
2 mi 





Figure 4.8 Discontinuity lattice configuration 


The flowgraph of the model shows an interaction of right- and left-propagating 
waves through the multi-section structure. We will characterize the lossy delay line 


for an m-th section by the transfer function 


(m) 
Din(z) = e7 [Ot © bm] 2-Tm (4.23) 
where al”) is the m-th section attenuation factor, l,, the section length, and Tj, is 


the discrete-time propagation delay per section. 


In summary, both equivalent circuit and ARMA/AR parametric models were 
presented for a cascaded microstrip step discontinuity. Losses were included in 
the circuit model by resistive -networks. Corresponding lossless reflection and 


transmission network functions were derived from the reactive networks. Next, 
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discrete transformations were performed on the network functions to derive the 
ARMA/AR digital transfer functions. Multiple microstrip step discontinuities were 
modeled by cascading the lossy equivalent circuit of a single step discontinuity. 
A lattice-like 1st-order ARMA model was introduced to describe a multi-section 
microstrip structure. In Chapter VI, transient analysis simulations of the equivalent 


circuit models will be compared with experimental measurements to validate 





Figure 4.9 Multi-section lattice model 


the accuracy of the circuit models. In the next chapter, we will develop several pa- 
rameter estimation algorithms, and introduce new system identification techniques 


for multi-section microstrip structures. 
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V. PARAMETER ESTIMATION ALGORITHMS 


The parameter estimation algorithms presented are based upon deterministic 
and stochastic methods. Figure 5.1 describes an unknown system to be identified. 
The observable input and output signals can be either deterministic or stochastic. 
The ARMA and/or AR model parameters of the unknown system are determined by 
the parameter estimation algorithms. Additionally, simulation results are presented 


to study the performance of the algorithms. 









SYSTEM TO BE 
IDENTIFIED 


IDENTIFICATION 


MEASURABLE ALGORITHM MEASURABLE 


OUTPUT 


SYSTEM MODEL 


Figure 5.1 Block diagram representation of the p’ .meter estimation problem 


The deterministic methods to be presented are the weighted least squares (WLS) 
algorithm, the network function approximation method, and the Schur algorithm. 
The WLS algorithm is a method for finding the filter coefficients of a p-th order 
recursive digital filter, -vhich gives an optimurn least squares approximation to a 
known sample impulse response. The network function method approximates an 
s-domain network function from the unit step response characteristics, such as rise 


and delay times, of a propagating pulse. Finally, the Schur algorithm estimates 
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the AR model parameters from a given sample impulse response. Our emphasis 
will focus on the deterministic methods because avaliable microwave measurement 
techniques have deterministic excitation sources. 

In the discussion of the stochastic methods we assume that a nearly white 
wideband microwave noise source excites the microstrip structure. When the noise 
source and the output random data are both available for signal processing, an 
ARMA parameter estimation algorithm based on a generalized Mullis-Roberts (M- 
R) criterion is used [Ref. 15]. However, if the input noise source cannot be measured, 
an estimate of the driving source can be obtained from the output process. This 
restriction in observable data precludes the direct application of the general estima- 
tors. Therefore, a suboptimal two-stage least squares ARMA estimation algorithm 
will be presented to solve this problem. 

In Chapter IV, we presented a lattice parametric model for the multi-section 
microstrip discontinuity. New algorithms will now be introduced that will estimate 
the ARMA/AR parameters for this model using reflected and transmitted time- 


domain measurements. 


A. DETERMINISTIC PARAMETER ESTIMATION 

1. Weigthed Least Squares (WLS) Algorithm 

a. WLS Algorithm Derivation 
Evans and Fischl [Ref. 43] presented a computational algorithm for es- 

timating the coefficients of a pth-order ARMA digital filter whose unit impulse 
respouse best approximates, in a minimum mean-square error sense, a prescribed 
impulse response of finite length. Here we derive a slightly modified version of their 
original least squares solution to estimate the ARMA parameters. 

The transfer function Harm a(z) of (3.2) can be expressed in terms of samples 


of the impulse response 


Harma(z) = ho + hyz7) + hoz? +... (5.1) 
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The first m samples of the estimated impulse response forms the m-vector 


ke) (5.2) 


ho 
h= _ (5.3) 
lisea 
The error vector is defined by 
e=h-h. (5.4) 


By long division of (3.2), we can relate the AR and MA parameters (a,b) to the 


| (5.5a) 


estimated impulse response h (a, 5) as 





where 
Po ag 
b=| | | anda=|° (5.56) 
ho 0 0 0 
[j= | "fo ae (5.5¢) 
hozae Rgacs ho 0 
By fips eh 
(H2] = host * a se , (5.5d) 
ot 
A least squares solution can be determined by minimizing 
a 1/2 
ll  (a,5) |] = bs ates? (5.6) 
1=0 
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with respect to the coefficients (a, }). 
Burrus and Parks [Ref. 44] have proposed a linear-equation error vector in the 


matrix formulation, given by 


d(a) = [H2] a (5.7a) 
where d 
i(a 
d(a) = at (5.78) 
dm—p(a) 


The error vector e(a, b) of (5.4) is related to d(a) by rearranging (5.7a) as 


d(a) = [B(a)]? h (5.8) 
where the m x (m — p) matrix 
ay 0 0 
Bia)=| ap a ... ap|- (5.9) 
0 ao «+. @Qy 
0 0... ay 


We can substitute & of (5.4) into (5.8), and recognizing from the lower partition of 
(5.5a) that 
[B(a)]? h=0, (5.10) 


we have 


d(a) = [B(a)]” e(a, 6). (5.11) 


An inverse relationship between the errors is formed using a transformation weight- 
ing matrix [W(a)]. Hence, the minimization problem of (5.6) can now be expressed 
as [Ref. 44] 

min |{ e(4, 6) || = min | [W(a)]d(a) | (5.12) 


where 


[W(a)] = [B(a)][B(a)? B(a)| (5.13) 














The computational algorithm begins with an initial estimate of 4). Equation 


(5.7a) can be partitioned as 


a(a)=[-f | [xX] @ (5.14) 
where 
hy-1 ho 
rn | 
PST Me se, (5.15) 
hm9 hm-—p-1 
—hy 
—h 
f=| (5.16) 
hm 
and the AR parameter vector 
a 
a 
a=| |. (5.17) 
Gp 


The least squares solution of a yields 
a) = [XTX] XT F (5.18) 
The initial estimate equally weights the error vector d(a). However, the error 


formulation can be weighted differently depending upon the initial estimate. Then, 


the weighted least squares solution at the i-th iteration is given by 


at) = [xtw (a) "w (a-)) x| es x 


|xtw (a) 'w (a) £| (5.19) 
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where the weighting matrix [w (a) is in terms of the previous estimate a!'~)), 


The iteration is continued until af) ~~ at), Using the optimal estimate for a, the 


weighted error e,,(a, 6) is calculated by 
€y(4,) = [w (a)] a(a) . (5.20) 


From (5.4), the estimate A(a, 6) is given by 


h(a, b) = h — ey(a, 6). (5.21) 


Finally, the numerator coefficient estimate vector 6 is calculated directly by 


— 


b= (H1] 


I> 


(5.22) 


b. Simulation Results 
The weighted least squares algorithm was implemented using a Fortran 
program. Several reference models were estimated beginning with an ARMA order 
(3,2), denoting AR order p = 3, and MA order g = 2. The transfer function is 


0.5 — 0.427! + 0.89277 


se 1 — 0.20z—1 — 0.252-? + 0.05273" 


(5.23) 
The actual reference model parameters and the WLS algorithm parameter estimates 
are shown in Table 1. All results are obtained after one iteration of the algorithm 


using 20 samples of the impulse response. 


Table 1. WLS SIMULATION (3,2) MODEL RESULTS 








We next consider a second reference model with order (4,2) having transfer function 


1+ 0.2z-! — 0.99272 


He) Sa a 
(=) 1 — 0.20z—-! + 0.62z~? — 0.1522-3 + 0.3016274 


(5.24) 


and the obtained results are given in Table 2. 


Table 2. WLS SIMULATION (4,2) MODEL RESULTS 





The final example is a strictly AR model of order (4,0) having a transfer function 


1 


(2) = ——— oOo 
(2) 1 — 0.20z-! + 0.6222-? — 0.1512-3 + 0.355274 


(5.25) 
The actual and estimated AR parameters are listed in Table 3. 


Table 3. WLS SIMULATION (4,0) MODEL RESULTS 





The above examples demonstrate the estimation accuracy of the WLS algorithm 


for (q¢ < p) order models. 


2. Network Function Approximation Method 
The following method, introduced by Elmore [Ref. 45], can approximate 
the normalized networ': function from the unit step response delay and rise times. 


Assume that the microstrip section is excited by a unit step and that its step 


51 











response v2(t), has been normalized as v2n(t). The microstrip group delay time is 


fi Vd_(t)dt 
oo = fe h,(t)dt (5.26) 


J sto : ‘ 


where v},,(t) = hp(t) is the normalized impulse response. The rise time is defined 


then defined as 


™= 


as 


9 


* 1/2 
2 / mina 
0 
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foe) 
fe =| fe 
0 


Figure 5.2 describes a network voltage output response resulting from a unit step 


input. The normalized network function is expressed as 
co 
= ‘| hy(t)e7 "dt (5.28) 
0 


—st 


Expanding e * in a power series in (5.28) yields 


0° 











As shown in Figure 5.3, the microstrip transmission line is described by 


M identical network sections iu cascade each with independent propagation delay 





Figure 5.2 A network N with its unit step response 


times Ty), Tg.,---, Tgaq> and rise times T,,,7;,, ...,T7; mu: ~hen the overall normalized 


network function is 


k=1 
M 2 2 
& T. 2 
=I out 5 (Get) - |. 
k=1 
M M OM 
3? T}, 2 
=1-s itt 5 On + Tok +25° Y= ToT 91 ee (5.30) 
=) : k=] l=k41 
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Finally. 
M ge M ee M 
H,(s)=1-s) tk + 5 pee (4 ae (5.31) 
i k=l * Lk= k=1 


where the total propagation group delay is 


M 
™%™= aor (5.32) 
k=1 
and the overall rise time for M sections is 
M 
P=) 7). (5.33) 
k=1 








Figure 5.3 Microstrip transmission line as M cascaded networks 


Assume that the propagating pulse is the response of some networ' to an unit step. 
Evidently, the rise time of the previous network(s) is T;;. Now the pulse continues 
through an unknown network which has a rise time 7,. The resulting pulse then has 


a total rise time 7-2. From (5.33), the expected rise time of the unknown network 
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is given by 








After estimating the network delay and rise times, the ARMA filter coeffi- 


cients of equation (3.11) are solved from the polynomial coefficients of (5.29), while 
the first-order AR filter coefficient is also related to 7, and T;, as follows: 


a = exp |-———,——"_ (5.35) 


where T denotes the sampling interval. 


3. Deterministic Schur Algorithm 
The autocorrelation function (ACF) of the normalized impulse response of 


equation (3.16) is computed as 
oO 
R(r) = i h(t)h(t+7)dt, for 7=0,1,2,...,N (5.36) 
0 


where 7 is the correlation lag. However, the sampled autocorrelation function will 


give a discrete correlation vector as: 


[R(0), R(1), R(2), ... ,R(N)] (5.37) 
Next, the ACF data is applied to the Schur algorithm. The Schur algorithm is a 
recursive algorithm that solves the AR model parameters (also called the partial 


reflection coefficients) [Ref. 14]. The algorithm begins by forming a generator 


matrix (G7) using the sampled autocorrelation vector as 


of = (70) FE) BR RN. (6.38) 


Shifting the first column down yields 


a= [5 88} ag 2 “fun?) oa 
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The Ist-order model reflection coefficient is computed as the ratio of the (2,2) and 
(2,1) terms of (5.39) as 


_ A) 


=——, 5.40 
A new generator matrix is then formed, given by 
GT = O(k,)GT (5.41) 
where 
1 1 —-k 
0(k,) = ———— |- u d (5.42) 
Similarly, the 2nd-order model reflection coefficient is 


* RO) — ki R(1)" 


The algorithm repeats until the desired order reflection coefficient is obtained. 
When the impulse response, h(t), describes a 1st-order AR transfer function, a sin- 
gle reflection coefficient will result. From this reflection coefficient the time constant 
of the impulse response is easily found. The computational efficiency is achieved by 
iterating the autocorrelation vector through an Nth-order AR lattice structure. In 


the next section, the Schur algorithm will be applied to an AR stochastic process. 


B. STOCHASTIC PARAMETER ESTIMATION 

The following algorithms are presented to estimate the ARMA/AR model pa- 
rameters when the observable input /output signals are stochastic processes. Assum- 
ing a white noise process as the driving input, the system to be identified produces 
an output process as shown in Figure 5.4. The efficiency of a particular estimation 
algorithm is governed by the choice of model selected and its order. In practice, we 
do not usually know a priori which model to choose. Once a model, either ARMA 
or AR, has been chosen, we must specify the model order. In choosing a model for 
the microstrip transmission line, we will select the ARMA/AR models proposed in 
Chapter III. 
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Autoregressive models are the most widely used models because the analysis 


algorithms for extracting the model parameters are found by solving a set of linear 


SYSTEM TO BE 
IDENTIFIED 


H(z) 
SYSTEM MODEL 





Figure 5.4 Stochastic system identification problem 


equations. When the AR modeling assumption is valid, these algorithms provide 
very good estimates of the model parameters [Ref. 33: p. 131]. Unfortunately, 
application of these algorithms to non-AR time series data usually results in poor 
quality estimates. Furthermore, if an AR estimator is applied to a process that is 
not AR, then the true AR model would be one of infinite order. Any finite order 
AR model will introduce bias errors from modeling inaccuracies. A trade-off must 
take place to choose order-p large enough to reduce the bias or to choose order-p 
small enough to reduce the estimation errors. 

In ARMA modeling, the best least squares estimate of the model parameters 
is, generally, a nonlinear function of the past observations. Nonlinear optimiza- 


tion techniques are usually computationally intensive and may not converge to the 
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global minimum. Two ARMA parameter estimation algorithms will be presented 


for stochastic data. 


1. Schur Algorithm - Revisited 
The AR filter coefficients (a),a2,a3, ...,@p) are related to the autocorrela- 


tion matrix of the data to be modeled and the noise variance (¢”) by the normal 


equations: 
Ryy(0) Ryy(-1)_....-- Ryy(—n +1) 1 a 
ee Fa) : . a ee _ = o . (5.44) 
Ryy(n — 1) “ee Ryy(1) Ryy(0) ap 0 


Several algorithms have provided efficient methods for solving the normal equations 
(Ref. 33: pp. 116-118]. The Levinson algorithm for solving Toeplitz systems of lin- 
ear equations is well-known in digital signal processing. However, the Levinson 
algorithm has been discovered to be computationally less efficient than an alter- 
native, the so-called Schur algorithm, because of a requirement to evaluate vector 
inner products (Ref. 12: pp. 6-22]. 

We have previously introduced the Schur algorithm for estimating AR model pa- 
rameters using deterministic signals. Assuming that the unknown system is driven 
by a white noise process, the AR model parameters can be computed from the 
output process using the Schur algorithm. The time average approximation of an 
autocorrelation function (called the sample autocorrelation) is defined as [Ref. 46: 


p. 27] 
N-1-k 


Ryy{n] = = Dy y[k+njy[k] for O<n<(N-1) (5.45) 
k=0 


where y[n] are tiie measured output signal samples of length N. A derivation of the 


Schur algorithm for stochastic signals is presented in Appendix A. 


2. Generalized Mullis-Roberts (M-R) Algorithm 
An ARMA parameter estimation algorithm based on a generalized Maullis- 


Roberts criterion has been proposed by Miyanaga et al. [Ref. 47: pp. 619-621]. 
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This algorithm is derived by forming forward and backward prediction models using 
available input noise and output data. Prediction errors are formed to satisfy specific 
orthogonality conditions which are similar to those found in AR modeling problems 
[Ref. 47: pp. 116-121]. Finally, autocorrelation and crosscorrelation functions 
are calculated and then the forward and backward ARMA filter coefficients are 
predicted from the correlation lags. A previous reference model is estimated using 
the M-R algorithm in order to compare parameter estimates obtained by the WLS 
algorithm. The transfer function of order (4,2), given by (5.24), was driven by a 


white noise sequence. Table 4 shows the ARMA parameter estimates obtained 


Table 4. M-R SIMULATION (4,2) MODEL RESULTS 


(WLS) Estimate | (M-R) Estimate 





from the M-R algorithm. 


3. A Two-Stage Least Squares Algorithm 

When the input noise is unobservable, the M-R algorithm cannot be applied 
directly. A modified two-stage least squares approach is now presented. The name 
two-stage least squares is due to the use of a large order AR least squares estima- 
tor followed by another ARMA parameter estimation algorithm. This approach 
is known to be suboptimal because the noise input is unobservable and must be 
estimated [Ref. 33: p. 320]. Nevertheless, we modify the two-stage least squares 
algorithm to estimate the specific ARMA model parameters of equation (3.11). The 
modification uses the M-R algorithm instead of the Levenison algorithm as sum- 


marized in Table 6. A reference model was estimated having a transfer function 
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of 
0.117227! 


1 — 1.451327) + 0.570z~? e-46) 


HarMa= 


at a sample period of 1.25 ps. Table 5 lists the results of this algorithm. The M-R 
algorithm (Step 3) produced very poor MA order coefficient results from the noise 


estimate. However, the AR order coefficients a; and az were used to calculate the 


Table 5. TWO-STAGE ALGORITHM RESULTS 


0.1172 0.1097 


-0.5700 -0.5701 


single MA order coefficient with favorable results. 






















C. MULTI-SECTION MODEL PARAMETER ESTIMATION 

In this section new algorithms are presented for the identification of the multi- 
section microstrip structure. The layer reflection coefficient estimation algorithm 
will approximate the number and location of cascaded microstrip step discontinuties 
from the measured reflection data. The layer-probing algorithm is an iterative 
impulse response estimation technique. The algorithm assumes that the multi- 


section microstrip structure can be described by the lattice AR parametric model 


of Chapter IV. 


1. Layer Reflection Coefficient Estimation Algorithm 
A single discontinuity between two lossy transmission lines having different 
values for the characteristic impedance Z,,, and propagation function 7 is shown in 
Figure 5.5. The m-th subscript denotes a single layer defined by equal propagation 
time delay. Assuming that a current wave iy is incident at the discontinuity from 
the m-th layer, we find that the incident current wave alone cannot satisfy the 


boundary conditions at the junction [Ref. 35]. 
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Table 6. SUMMARY OF TWO-STAGE LEAST SQUARES ALGORITHM 


Given the output process, y[n], of length N, the sample interval, T, 
and L > N 


Step 1: Predict the Lth-order AR parameters (4;) 


Using the Levinson algorithm [Ref. 33: pp. 213-214] with output, (y[n]) 


Step 2: Estimate the noise input (é[n]) 


Estimate AR filter coefficients (a,) and (a2) 
Use the Mullis-Roberts criterion algorithm [Ref. 15] 


where (é[n]) => input and (y[n]) > output 


Solve (Cw T) 
Substitute (a2) of Step 3 => (3.11c) 


Solve (wo\/1 — ¢? T) 
Substitute (a) from Step 3 = (3.11b) 
Substitute (Cwo T) from Step 4 > (3.11b) 


Solve MA filter coefficient (6;) 
Use (wo T) from Step 4, and 


(wo./1 — ¢? T) from Step 5 to solve ¢ 
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Hence, a reflected current wave ie” is generated. Using the boundary conditions 


at the impedance mismatch, we then write 


tig = Sa +i), 


where 


ity = i) em lm 


oP fo) rm bo, 


d=0 } d=C 
| gle) gC) pf+) 
m m m+1 
SS 
| | 
es ue Z ntl Mime 


Figure 5.5 Impedance discontinuity equivalent circuit 


The voltage and current wave relationships are defined as follows: 


(4) tae 
tm = SoS" 
m 
(-)) _ of” 
= 7 
BCA) ae Dior 
m+1 Zm+1 
Pe eas Ut 
me Zm+1 


(5.47a) 


(5.476) 
(5.47c) 


(5.48a) 
(5.486) 
(5.48c) 


(5.48) 





By substituting (5.48) into (5.47) we have 


ia Ba “ 
vt” — of” = (oft), - of ),). (5.49) 


The instantaneous voltage at the mismatch is defined as 


(-)! 


ota) = ult) — Um 
= oft) - ie (5.50) 


Taking (5.49) and (5.50), we solve for the incident voltage wave as 


1 ~ m 
ve) =o), (+) +00), (2) (5.51a) 


T 


where the reflection coefficient p,,, and transmission coefficient 7, at the mismatch 


are 


Zm+1 < Zim 
== 5.516 
and 
ps 22m 
ne Zm+) zs Zm 
Tm = 1+ pm- (5.51c) 
While the reflected voltage wave at the m-th layer is 
(-y _ 4) fPm\ -) fb 
oath (2) () 
= vn” (om) — vy (1 — pm) (5.52) 
Using (5.51a), the voltage wave incident at the (m + 1)-th layer is 
vt) = uff" (1 + pm) — 0), (pm): (5.53) 
Combining (5.52) and (5.53) into a scattering matrix yields 
+ 
cee ie om (5.54) 
ub 4 Pm 1 Pm Un+4l 
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Next, we normalize the voltage waves with respect to their characteristic impedance 


as 
(+)! 1 (+) 
Vv —_ _* |Um 2 
head ee keg ee 
and 


sc) i> gt) 
Cy | = eS |.0'I- (5.555) 
Vin+1 EL LY, 


m+1 


Using (5.54) and (5.55), the normalized forward recursion matrix yields 


+ 
[tes | = ee E =] [78 . (5.56 
vO | Via lem 1 YL 


The normalized incident and reflected voltages in (5.56) can be related to the prop- 


agation function (7m) as 


—(@m+j8m)lm (+) : 
i s som tienes cy . (5.57) 


Therefore, (5.56) can be expanded as 


—Bm ln (+) 
x E ebm 2 Ral (5.58) 


For a lossy dispersive microstrip transmission line the attenuation am, and phase 


Bm terms are functions of frequency. 
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Hence, the total attenuation per m-th layer is described by 


Om(f) = ae(f) + aa(f) (5.59) 


where the conductor loss a,, and the dielectric loss ag were given by (2.8b) and 
(2.10), respectively. Whereas, the phase term is related to the frequency-dependent 


effective microstrip permittivity es ¢(f) as 


Rafe ee) erie (5.60) 


Both the attenuation and phase terms are approximated as constants at the m-th 
layer. Furthermore, the phase term {,, is proportional to the m-th layer propagation 
delay time [Ref. 35] as 


lin 


Ugm 
dBm La 
= Law (@) (5.61) 


where vg,, is the group velocity in the m-th layer. 

Next we develop the discrete-time recursion equations that solve the inverse 
scattering problem. Taking the z-transform of the recursion matrix (5.58) and using 
(5.61) yields 








Vangie) -| v2) ie 
Vere: - a Vin (2) ee) 


where ; 
Ae 2) — pm || e~amlm 
[9] = Fel Spat || Fee 


The variable z~' represents the two-way propagation time delay of the m-th section, 
: 1 ; eee A 

while z~? represents the one-way time delay. Using (5.58), the (m + 1)-th laver 

reflection coefficient is computed by 


a 
Pm4i = aes (5.63) 


m+1 
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A layer characteristic impedance profile (in ohms) is determined from the layer 


reflection coefficient sequence as 


(5.64) 


1 m 
ya ae ae 


1 — pm+i 


where the initial impedance Zo is the measured input characteristic impedance. In 
order to compensate for losses, the algorithm must have prior knowledge of the 
total attenuation per layer a. Unfortunately, this loss information is not available 
for the inverse scattering problem, but the overall attenuation per structure can be 
estimated from the measured incident, reflected, and transmitted signals. Under 
this assumption, we can obtain an average attenuation factor for each layer. A 
method to estimate the attenuation factor from measured results will be introduced 


in the last section. 


2. Layer-Probing Algorithm 

The multi-section microstrip structure is modeled as cascaded linear time- 
invariant systems. Figure 5.6 describes the first iteration of the layer-probing algo- 
rithm. The algorithm begins by exciting the structure with an impulsive signal v.[n]. 
A minimum-phase transfer function 7\(z) is estimated for the first discontinuity us- 
ing the WLS algorithm. Next, an inverse filter A;(z) of the transfer function 7}(z) is 
implemented and excited by the original impulsive signal. The inverse filter output 
now becomes a new excitation signal. This new input is convolved with the first 
discontinuity transmission impulse response, resulting in an impulsive signal vp[n] 
which probes the second discontinuity. The subsequent reflected wave v,{n] now 
contains the reflected impulse response of the second discontinuity. We can ignore 
that portion of the reflection response caused by the first discontinuity, and record 
the subsequent data. The recording window is defined by a time greater than the 
two-way propagation delay time to the first discontinuity. Next, the windowed re- 
flection response is deconvolved across the first discontinuity transmission impulse 
response in order to observe the second discontinuity reflected impulse response. 


Similarly, another AR model is estimated from this reflected impulse response, and 
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then we repeat the procedure. Practically, the excitation signal is bandlimited by 


the instrumentation, and this prevents the input from being an ideal impulse. 


f 


: ne ie Cel, f=aL 
s | 
| 


Inverse Filter 


Vv 
i 
| 


| 
Vo {n] = 


Ist Discontinuity 2nd Discontinuity 


Figure 5.6 Layer-probing of a unknown discontinuity 


Also, the microstrip dispersion causes pulse spreading. As a consequence, when two 
discontinuities are separated by less than the time-dispersive spread for the propa- 
gating pulse, multiple reflections may introduce destructive interference. Therefore, 
when the reflected impulse response is significantly distorted by destructive inter- 
ference, the inverse scattering problem cannot be solved. 

We will attempt to overcome this problem by measuring the output trans- 
mission response resulting from the last impulsive probing signal. The remaining 
unknown discontinuities will be described by an M-th order AR model. A technique 
is required for estimating an appropriate model order. It is generally observed that 
M<€N where N is the data segment length. The best order is typically within the 
range from N/20 to N/5, but for particularly short data segments the best order 
can be as large as N/3 or even N/2 [Ref. 48]. Consequently, we do have a priori 
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knowledge about the order by assuming the model order and the detected num- 
ber of discontinuities are equal. As discussed earlier, the layer reflection coefficient 


algorithm can estimate the number of discontinuities within the unknown structure. 


The transmission response now reveals sufficient transient information to 
model the remaining unknown sections. Using the WLS algorithm, the transmission 
impulse response is modeled by an ARMA transfer function T(z). Suppose for 
simplicity, and not too much loss of generality, that we let the MA model (numerator 
term) become a gain constait, so that 

I(z)= ra (5.65) 


and the transfer function can be rewritten as 


M 
T(z) = |] T(z) (5.66) 


r=1 
where 
pil 
T;(z) = 
r(2) 1—@, 27} 


and a; is the AR model parameter. As can be seen, the individual T,(=) functions 
are lst-order AR transfer functions which model each discontinuity. Simulation 
results are performed in Chapter VI to demo \strate the performance of the layer- 


probing algorithm. 


3. Model Attenuation Factor Estimation 
The total attenuation factor can be estimated by considering the energy ]«ss 
of the microstrip structure. An incident transient pulse delivers a finite amount 
of energy to the microstrip transmission line. Discontinuities in multi-section mi- 
crostrip structure will introduce multiple energy reflections. We can only measure 


the refle :t 1 and transmitted voltage waveforms with respect to the input/output 
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measurement ports. Thus, we will treat the lossy multi-section structure as a sin- 
gle cascaded microstrip step discontinuity. Referring to Figure 5.7, the cumulative 
energy loss is given by 


Exoss = Ey) — Ey) — Ey” (5.67) 
where 


EM) is the cumulative reflected energy measured at the input port (1), 

EB is the cumulative transmitted energy measured at the output port (2), 
and 

EO) is the cumulative incident energy measured at the input port (1). 


The cumulative energy values are calculated by [Ref. 31: p. 34] 
N 
E =) |v[n])? (5.68) 
n=0 


where v[n] is the sampled voltage measured at the respective input/output ports. 
Next, the energy loss is described in terms of an imaginary discontinuity boundary 
(D) as 
E,oss = (Ey - EY) + (ED _ Et”) + (Ey - E}) (5.69) 
where 
EW) _ EM e201 h 
EY = ET) 2a h 


ja = EX) gar Ip 


Assuming a symmetrical microstrip geometry with respect to the boundary. we will 
assume Q] = a9 and 1, = ly. Equating (5.67) and (5.69) yields the half-symmetrical 


loss ( in nepers ) as 


(r) (t) 
ay i= an Bice . (Np) (5.70) 
(1) 
4| Ey 
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The half-symmetrical attenuation factor will be used to design the resistive 7z- 


networks in Figure 4.5. 





(r) | - 


<—_ <= 
a ae ear ae 


Figure 5.7 Cascaded step energy flow model 


This chapter presented the use of several ARMA/AR model parameter esti- 
mation algorithms for the characterization and transient response modeling of mi- 
crostrip transmission lines. Both deterministic and stochastic approaches were dis- 
cussed. Newly developed algorithms to improve the identification of lossy microstrip 
lines consisting of multiple discontinuities were also derived. These algorithms were 
the layer reflection coefficient estimation algorithm, and the layer-probing algorithm. 
The purpose of these algorithms is to estimate the appropriate ARMA/AR model 
parameters from available tiine-domain measurements. In Chapter VI, the algo- 
rithms’ performance will be studied using simulated and experimentally measured 


data. 
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VI. EXPERIMENTAL AND SIMULATION RESULTS 


The purpose of this chapter is to verify the proposed microstrip models using 
both experimental and simulation results. The lumped-distributed equivalent cir- 
cuit and parametric models presented are now simulated for an integrated circuit 
(IC) interconnection, a cascaded microstrip step discontinuity, and a multi-section 
microstrip step discontinuity. The microstrip models are verified by comparing pi- 
cosecond time-domain measurements with transient analysis and parametric model 
simulation results. Finally, the performance of the ARMA/AR parameter estime- 


tion algorithms is studied using simulated and measured data. 


A. INTEGRATED CIRCUIT INTERCONNECTION RESULTS 


1. Optoelectronic Measurement Techniques 

D.H Auston and A.M. Johnson of Bell Laboratories first demonstrated the 
generation of picosecond pulses on silicon substrates in 1975 by using optically trig- 
gered gating [Ref. 49,50]. They used quasimetallic photoconductivity produced by 
the absorption of short optical pulses in silicon to both turn on and off a switching 
gap. Figure 6.1 shows an optoelectronic switch using this technique. The microstrip 
transmission line is fabricated on high resistivity (p = 104Q — cm) silicon and con- 
tains a gap which forms the switch. Optical pulses generating mobile carriers within 
the silicon produce the switching action. The work on silicon has been extended by 
W.R. Eisenstadt and R.B. Hammond to produce picosecond optoelectronic switches 
on bulk silicon [Ref. 51]. Their effort was focused towards the production of high 
speed photoconductors which can be integrated along side other devices on an in- 


tegrated circuit using standard VLSI fabrication techniques. 


The Photoconductor Circuit Element (PCE) is the best laser driven optoelectronic 
transducer for high-speed sampling applications to integrate on a silicon substrate 
[Ref. 52]. The primary goal of the PCE is to achieve an ideal impulsive sampling 
gate. 


Optical Pulse 






Microstrip Line 


Ground Plane 


Figure 6.1 Diagram of D.H. Auston’s picosecond photoconductor 


2. Picosecond Pulse Propagation Measurements 

D. Bowman [Ref. 53: pp. 82-94] performed a series of picosecond time- 
domain measurements using a polysilicon PCE test structure. The structure shown 
in Figure 6.2 illustrates the IC fabrication of the PCE sampler on the microstrip 
transmission line. The length of microstrip beyond the last PCE sampler prevents 
any reflections from the open stub end from returning in time to interfere with 
the measurement. The experimental sampling system uses a colliding pulse mode- 
locked (CPM) ring dye laser with a nominal 200 femtosecond optical pulse to excite 
a 3.5 GHz bandwidth eiectrical pulse. The half power bandwidth of the electrical 
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electrical pulse is determined from the magnitude spectrum using the Fast Fourier 


transform. Table 7 lists the physical properties of the microstrip test structure. 


PCE SAMPLER 









NS NN 
IN ! 1 1 Microstrip Line + 1 


4 — 
IN-LINE PULSING PCE 





Figure 6.2 Polysilicon PCE microstrip test structure 


a. Sampled Waveform Measurements 
A transient pulse was measured at different distances along a 100 ym 
wide microstrip at a sampling interval of 1.25 ps. Figure 6.3 shows the measured 
waveforms with magnitudes normalized to the peak of the first sampler. Successive 
curves are the sampled mvasurements at the second (500 ym), third (1000 zm), and 
fourth (2000 um) PCE sampler. 
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Table 7. IC MICROSTRIP PROPERTIES 


Material Properties 


Substrate Resistivity (Si) (ps) 
Conductor Resistivity (Al) (p-) 
Conductor Length ()) 
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Figure 6.3 High frequency waveform propagation measurements 
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The measured responses show the expected dispersive propagation with high losses. 
This is evident from the increasing rise times of the pulses (pulse spreading) and 


smaller peak magnitudes with increasing propagation distance [Ref. 53: pp: 84-85]. 


b. Propagation Velocity 

The rate at which a pulse propagates down an IC interconnection is a very 
important parameter, particularly for circuits involving critical timing applications 
or pulse synchronization. At the high frequencies present with picosecond pulses, 
propagation velocity varies with frequency and affects the shape of the pulse. 

The relative time at which the peak magnitude of the pulse occurs can be 
used to determine the velocity. The group velocity and not the phase velocity of the 
pulse is actually measured by observing the relative delays of the peak magnitude. 
Figure 6.4 reports the results of the measured propagation delays at the various 
time differences between the peak magnitudes for each of the three samplers. The 


theoretical propagation delay is calculated by [Ref 18: p. 62] 


€ 1 
ea ma (6.1) 





for the different lengths of the microstrip. What is interesting about the results 
is the almost identical slope for both the measured and theoretical delay. This 


indicates a relatively constant velocity of ~ 112-114 ym/ps. 


c. Propagation Loss 

The decrease in the peak magnitude of the propagating pulse with distance 
seen in Figure 6.3 is evidence of the significant losses present on silicon substrates. 
Figure 6.5 shows the loss in decibels as a function of propagation distance. The 
theoretical loss combines the numerical results for conductor loss (2.8b) and dielec- 
tric loss (2.10) at three microstrip distances. The measured propagation losses were 
obtained from relative cumulative energy results. The theoretical values reveal a 
loss of 2.4 dB/mm while the measurements indicate a 3.6 dB/mm loss is actually 


occurring. A possible explanation for the discrepancy between the theoretical loss 











and the measured loss is that (2.8b) yields a somewhat low result because surface 
roughness is not considered, and the effective microstrip permittivity value used in 


(2.10) actually varies with frequency [Ref. 18: p.90]. 
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Figure 6.4 Pulse propagation delay on IC interconnection 


3. Equivalent Circuit Simulations 
The lumped-parameter and lumped-distributed equivalent circuits, derived 
in Chapter II, are developed for a 500 ym length of microstrip in the PCE test 
structure. Numerical results for the circuit model are presented in Table 8. Several 
comments can be made about these theoretical calculations. The static effective mi- 
crostrip permittivity significantly influences all subsequent calculations. It assumes 


a TEM mode of propagation on the microstrip which is frequency independent. 
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In the equivalent circuit model the lumped-elements are evaluated at a maximum 


frequency of 3.5 GHz which corresponds to the bandwidth of the propagating signal. 
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Figure 6.5 Pulse propagation loss on IC interconnection 


At this frequency the physical microstrip length should be less than 4 mm (* A, /8). 
Consequently, the equivalent circuit models may lose their accuracy if significant 


higher order modes are excited by microstrip discontinuities. Figure 6.6 shows 
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Table 8. gh a atrodad CIRCUIT THEORETICAL RESULTS 


Circuit Parameter 


ae 
oa) 
Shunt Capacitance < (2. 18a) 
a 
ORE 


Total Loss (az) a 1.213 dB 


the schematic diagrams of the equivalent circuit, models. 

















74.23Q 10.42 0.332 nH 


74.239 





74.23Q 


(a) 
Figure 6.6 (a) Lumped-element equivalent circuit 


In modeling the IC interconnection using equivalent circuits, it is assumed that 


impedance matching networks are used to prevent reflections. In Figure 6.6(a) a 
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propagation loss of 1.213 dB is modeled by a resistive r-network and a propagation 


delay of 4.46 ps is described by the reactive LC network. Transient analysis simu- 


10.40 


ame 4.46 ps 


74.232 





(b) 
Figure 6.6(b) Lumped-distributed equivalent circuit 


lations of the circuit models are performed using the PSPICE electrical circuit simu- 
lation program (Ref. 17]. The PSPICE network listings for both equivalent circuits 
are given in Appendix B. A piecewise-linear approximation of the normalized refer- 
ence pulse, measured at the first PCE sampler, acts as the PSPICE voltage source. 
Waveform comparisons between the second PCE sampler measurement and the 
circuit model simulaticns are shown in Figure 6.7. The general shape of the simu- 
lations are representative of the measured data; however, there are some significant 
discrepancies. Two possible causes are suggested. First, on the decaying transient 
the measurement contains a second peak caused by a reflected signal from the back- 
plane of the substrate. This phenomena is not modeled by the circuit model because 
of its complex wave nature [Ref. 53: p.90]. Second. the equivalent T-network circuit 
model simulation indicates a decreased peak magnitude. This observation can be 


explained as energy lost by a reflected voltage caused by the shunt capacitance. 











Next, we derive an equivalent circuit model for the 1000um length of microstrip. 
Using the same reactive 7-network of Figure 6.6(a), we calculate the new shunt 


capacitances and series inductance as 0.06 pF and 0.663 nH, respectively. The 
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Figure 6.7 Waveform simulations of the equivalent circuit models 


measured propagation loss of 3.891 dB is used to evaluate the resistive -network. 
Figure 6.8(a) compares the PSPICE transient analysis simulation with the measured 
response. The simulated results tend to underestimate the rise time of the measured 


pulse. 
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This error may be attributed to the breakdown in the lumped-element equivalent 
circuit assumption as the physical length of the microstrip increases. Therefore, we 
can improve the model by cascading two reactive x-network sections. A PSPICE 


simulation was performed using the higher order circuit model. 
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Figure 6.8(a) Single LC network simulation at the 1000u4m PCE sampler 


Figure 6.8(b) shows a significant improvement in the waveform agreement, especially 
in the rise times. However, as the number of reactive lumped-elements increase, the 


model order also increases. 
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Then the mathematical relationships among the model parameters and the mi- 


crostrip physical properties will become more complex. 
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Figure 6.8(b) Cascaded LC network simulation at the 1000um PCE sampler 








4. Parametric Model Simulations 
The following ARMA digital transfer function is computed for the 500 wm 
microstrip test section using (3.10) and (3.11) at a sampling interval of 1.25 ps: 


0.1172 27} 


fn Es 6.2 
1 — 1.4513 2-1 + 0.57 27? (6.2) 


Harma(z) = 


The normalized reference pulse, measured at the first PCE sampler, is attenuated 
by 1.213 dB and then iterated through a difference equation of (6.2). Similarly, an 
AR digital transfer function is evaluated using (3.19) as 


0.4295 
Har(z) = 


1 — 0.5705 2-1)" 63) 


The normalized reference pulse is again attenuated by 1.213 dB and iterated through 
a difference equation of (6.3). Figure 6.9 shows the ARMA and AR model simulation 


results as compared with the measured pulse. 
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Figure 6.9(a) 1st-order ARMA parametric model simulation 
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Simulated pulses are shown by solid lines while the measurements are indicated 
with dash lines. Both simulated waveforms have properly accounted for the loss 
mechanisms on the 500m microstrip section; however, the 2nd-order ARMA model 
result shows an outstanding waveform agreement as compared with the lst-order 
AR model. In the next section we will focus on the model parameter estimation 


algorithms using both simulated and measured data. 
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Figure 6.9(b) 1lst-order AR parametric model simulation 


5. Model Parameter Estimation Results 
6 WLS Algonthmn: 
The weighted least squares (WLS) algorithm of Chapter V can be used to 
estimate model parameters from a given sample impulse response. First, we must 
obtain an accurate impulse response from the available input/output measurements. 


The impulse response was obtained by deconvolving the output data, measured at 
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the second PCE sampler, with the measured normalized input data. The original 


measured data, sampled at 1.25 ps, was decimated to a sampling interval of 5 






ps to improve the frequency resolution of the FFT during deconvolution. Next, 






the resulting sample impulse response was applied to the WLS algorithm. The 






estimated ARMA model parameters are listed in Table 9. Figure 6.10 compares the 






Table 9. ARMA PARTE USING Wie ALGORITHM 
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estimated impulse response to the measured and theoretical impulse responses at a 






5 ps sampling interval. 
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Figure 6.10 Comparison of ARMA model sample impulse responses 











b. Stochastic Methods 
ARMA and AR sample impulse responses were simulated using a PSPICE 
transient analysis of the equivalent circuit models at a sampling interval of 1.25 ps. 


Stochastic output processes are generated using the convolution summation 


N 


y[n] = So ofk]h[n — &] (6.4) 


k=0 


where v[n] is an input Gaussian white noise sequence and h(n] is the simulated 
sample impulse response. The white noise sequence was obtained using a pseudo- 
random number generator program. The output process y[n], having a record length 
of 1800 samples, and the original white noise sequence were applied to the M-R 
algorithm. The ARMA digital filter coefficients were estimated for a (2,1) order 
model. Table 10 compares these estimates to the theoretical parameter values at a 


1.25 ps sampling interval. 


Table 10. ARMA PARAMETERS USING M-R ALGORITHM 


Theoretical Estimated 





The estimated results compare favorably with the theoretical model parameter val- 
ues. There is, however, a limitation in using the stochastic approach. First, the 
M-R algorithm requires input/output data to be simultaneously sampled to provide 
accurate ARMA parameter estimates. Second, microwave measurement techniques 
are usually based on frequency-domain scattering parameters, and high-speed pi- 
cosecond sampling techniques are not commercially available. 

The Schur algorithm is used to estimate a lst-order AR model. Initially, a 
lst-order AR sample impulse response was simulated using PSPICE at a sampling 


interval of 1.25 ps. Using equation (5.45), a sample autocorrelation sequence is 
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computed from the impulse response data. Then the Schur algorithm is used to 


solve the AR reflection coefficients, and the normalized AR transfer function yields 


1 = 
Hapr(z) = To0s71* (6.5) 


Since the impulse response was simulated using a 1st-order AR model a single Schur 


reflection coefficient should be expected. 


c. Elmore Method 
Network functions are now estimated from rise and delay time measure- 
ments. Figure 6.11 describes two consecutive PCE sampler voltage measurements 
denoted as v2(t) and v3(t). A two-port network describes the 500um length of 


microstrip. 





Figure 6.11 Cascaded network functions 


The Elmore rise time (7;) is defined as the reciprocal of the slope of the tangent 
drawn to the response curve at its half-magnitude point [Ref. 46]. The delay times 
are estimated at consecutive peak or half-magnitude points. Referring to Figure 
6.12, a delay time of 4.41 ps is obtained by averaging the peak magnitude delay 
(4.48 ps) and the half-magnitude (4.33 ps) delay times. A rise time of 2.06 ps is 
computed using (5.34) for T;2 = 7.73 ps and T;3 = 8.00 ps. 
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A 2nd-order network function is evaluated as 


99.3836 x 102! 


H =e 
(s) = 37935080 x 108s + 99.3836 x 102! 


(6.6) 


The characteristic impedance of the microstrip must be determined prior to solv- 
ing for wo and ¢ in (3.5). The estimation of the characteristic impedance will be 
addressed in the next section. Finally, the ARMA model parameters are calculated 
at the sampling interval of 5 ps using equations (3.10) and (3.11). 
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Figure 6.12 Elmore delay and rise time graphical estimates 


Equation (5.35) solves the 1st-order AR filter coefficient in terms of the delay and 
rise times. Both the estimated and theoretical ARMA/AR model parameters are 
compared in Table 11 at a 1.25 ps sampling interval. 


88 














The Elmore estimates for the AR model parameters, including the denominator AR 
filter coefficients of the ARMA model, all show a slight increase over their theoretical 


Table 11. MODEL PARAMETERS USING ELMORE 


ARMA Parameter 


[ARPeranctee [| 










1, 









values. This result will explain the increase in effective microstrip permittivity 


estimates presented in the next section. 


d. Estimation of Physical Microstrip Properties 
The characteristic impedance at each PCE sampler is con:puted using 


(2.29) and (2.30). Table 12 lists the calculated characteristic impedance estimated 





at each PCE sampler. The average impedance of 74.32 2 compares very favorably 


to the actual microstrip characteristic impedance of 74.23 2. 


Table 12. CHARACTERSITIC IMPEDANCE ESTIMATES 


The effective microstrip permittivity (€esf) is directly calculated from the 















model parameter estimates using (3.22) and (3.25). Several theoretical expressions 
have been cited in the literature for the effective microstrip permittivity. Here, we 
have computed a value of 7.082 using equation (2.2). However, other closed-form 


expressions will produce different results. For example, another expression for the 
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effective microstrip permittivity, given the characteristic impedance, is [Ref. 18: 


p.44] 














where Zo = 74.23 M and e, = 11.7 for silicon substrate. Table 13 lists the effective 
microstrip permittivity values that were computed using the ARMA/AR model 


parameters obtained from the different estimation algorithms. 


Table 13. EFFECTIVE PERMITTIVITY ESTIMATES 


aR 


These results indicate that the best model parameter estimators are the WLS and 
















Schur algorithms. These algorithms provide optimal estimates in a least squares 
sense, while the Elmore method was extremely dependent upon obtaining accurate 
rise and delay time measurements. Precautions should be made whenever we com- 
pare the estimated results to the theoretical values. For example, the accuracy of 
the theoretical expressions is dependent upon the shape ratio w/h range. In all 
cases the shape ratio wil] be accurate to +1 percent. For narrow lines (w/h < 1.3), 
the effective microstrip permittivity has an error range +0.5-0.0 percent. When 
calculated using (6.7), eer is accurate to +1 percent [Ref. 18: pp. 45-46]. Ad- 
ditionally, the effective microstrip permittivity is in fact frequency dependent, and 


the values tend to be slightly higher than those given by (2.2) and (6.7). 
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B. CASCADED MICROSTRIP STEP DISCONTINUITY RESULTS 

The equivalent circuit and parametric models of Chapter IV are validated by 
comparing model simulations with experimental measurements. Figure 6.13 de- 
scribes a cascaded microstrip step discontinuity structure that was fabricated on 
G-10 epoxy dielectric material. The physical dimensions of the cascaded step dis- 
continuity were restricted by the available photolithographic equipment and the 
dielectric substrate. A half power bandwidth of 400 MHz for the measured refer- 
ence pulse was selected as the maximum design frequency. This frequency is used in 
the theoretical calculations. Table 14 summarizes the theoretical microstrip results 


obtained from the closed-form expressions. 





Figure 6.13 Cascaded microstrip step discontinuity 


1. Validation of Equivalent Circuit Models 
The lumped-element values in Table 14 are associated with the equivalent 
circuit models of Figures 4.3(c) and 4.4(b). A transient analysis simulation, using 


PSPICE, is performed for each equivalent circuit. A piecewise-linear approximation 








of the measured input pulse was used in the transient analysis simulations. A 
half-symmetric attenuation factor of 0.0348 nepers is computed from the measured 
incident, reflected, and transmitted cumulative energy. Separate resistive attenuator 
x-networks are designed to model each half-symmetric loss. The shunt resistors 


(2.316 kQ) are given by [Ref. 37: pp.186-192]: 








N+1 
Ry = Zo (74) (6.8a) 
and the series resistor (1.41 ]) by 
N? —| 
R,=Z 6.8) 
»= (SS) (6.86) 


where the loss in decibels is 20log;) MN. Transmission and reflection output. re- 
sponses are simulated. These simulations are compared to the measured responses 


in Figures 6.14(a) and 6.14(b), respectively. The measured transmission 
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Figure 6.14(a) PSPICE simulated transmission response 
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and reflection responses were obtained using the inverse FFT of the measured scat- 


tering ~..rameters. 


2. Parametric Model Simulations 
Using equations (4.14), and (4.16) — (4.18), yields a 3rd-order ARMA trans- 
fer function 


_ 6.776 x 10° — 7.295 x 10°27} + 4.323 x 10°27? 


T3(z) = 1 — 1.702=-! + 1.59642-2 — 0.762:-3 (6.9) 


at a sampling interval of 40 ps. A sample impulse rcsponse is generated using the 
difference equation of (6.9). The sampled incident pulse is then convolved with the 


sample impulse response data to produce an output transmission response. 


Table 14. CASCADED STEP THEORETICAL RESULTS 
~ Mires Paras 

"Wide Width (em) Segment | 
$7.00 


[Rerrow Width (0:97 em) Segment [| 
229) 
Series Inductance (Ls) 


Step Discontinuity Circuit Parameters eee 
Shunt Capacitance (Ca) 














Microstrip Parameter 
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The normalized transmission response is compared to the measured output in Figure 


6.15. Next, the Ist-order ARMA reflection and transmission transfer functions are 
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Figure 6.14(b) PSPICE simulated reflection response 


obtained using (4.20b) and (4.20c), respectively, as 


_ 0.746 (27! - 1) 


(2) = ————_———— 6.10 

MMe) = 9493271 aa 
and ’ 1 

0.254 (1 + 27 

15s Racawns 6.100 

NM) = 0 oie ae 
at a sampling interval of 200 ps. 
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Figures 6.16(a) and 6.16(b) compare the ARMA reflection and transmission re- 


spouse simulations with the measured responses. The results show both the lossless 
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Figure 6.15 3rd-order ARMA model transmission response 


and lossy (0.0348 nepers) simulated output waveforms. When the attenuation loss 
is included, the waveform agreement improves, especially with respect to the peak 
amplitudes. 

Finally, a 1st-order AR model is given by (4.21) as 


0.68 


Tart?) = 79 507s-1 


(6.11) 


at a sampling interval of 200 ps. The normalized AR transmission response is 


compared to the measured response in Figure 6.17. 
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Figure 6.16 1st-order ARMA (a) reflection response and (b) transmission response 

















3. Estimation Results 
a. Estimation of Reflection Impulse Response 
An appropriate parametric model for the cascaded microstrip step dis- 
continuity is estimated from only the measured incident and reflected responses. 
An estimate of the reflection sample impulse response is obtained by deconvolving 
the measured reflection response with the incident pulse. Figure 6.18 shows the 


estimated reflection sample impulse response at a sampling interval of 200 ps. 
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Figure 6.17 1st-order AR transmission response 


b. Estimation of Model Parameters 
The WLS algorithm estimates the 1st-order AR model parameter of equa- 
tion (4.21) using the previously estimated reflection impulse response. The time 
constant (C,Zo2) is solved as 291.51 ps using an estimated 1st-order AR filter co- 
efficient of 0.5035. 
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Next the lst-order ARMA model parameters are calculated by substituting this 
time constant value into equation (4.20c). Therefore, the unknown transmission 
response can also be approximated from observing only the incident and reflected 


data. Finally, the ARMA and AR model parameter estimates are given in 
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Figure 6.18 Estimated reflection sample impulse response 


Table 15. This approach will be exploited by the layer-probing algorithm in the 
next section for the system identification of a multi-section microstrip structure 


(Ref. 54]. 





Table 15. MODEL PARAMETER ESTIMATES 


ARMA Parameter 


[aRPaancer [id 


es 
















C. MULTI-SECTION STEP DISCONTINUITY RESULTS 

As we saw earlier, a cascaded step microstrip discontinuity can be modeled by 
an equivalent circuit or by an AR parametric model. We will now extend these 
models for a multi-section microstrip step discontinuity. 

1. Equivalent Lumped-Distributed Circuit Model 

a. Theoretical Results 
The capacitive shunt equivalent circuit model, as shown in Figure 4.4, 

will be used to characterize each cascaded step discontinuity. The high impedance 
(narrow microstrip width) lines which connect each discontinuity are modeled as 
lossless distributed transmission-line sections using (2.24). The propagation delay 
time associated with each transmission line is determined from the physical length 
and the frequency-dependent effective microstrip permittivity. Figure 6.19 describes 
the multi-section microstrip structure which will be modeled and simulated using 
PSPICE. Three different line width dimensions are referred to by numbered line 
segments. Specific equivalent circuit numerical results are listed in Table 16 for 
each line segment. The theoretical calculations use a maximum design frequency 
of 4 GHz and a substrate relative dielectric constant ¢, of 4.3. Frequency disper- 
sion is incorporated into the circuit model by using frequency dependent effective 


microstrip permittivity €e¢7(f) and characteristic impedance Zo(f) values. 
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The characteristic impedance is calculated using (Ref. 18: pp. 81-82] 
376.6 A 
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Oe Gen) ve 


where 
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Figure 6.19 Multi-section microstrip discontinuity test structure 


The theoretical propagation losses are given in Table 17 for each line segment. The 
total segment loss in dB/A, is the sum of each segment dielectric ag and conductor 
a, loss, given by 

_ 27.3 €, (cess - 1) tan é 


2.10 
wea ay) 

_ 0.072VF Ay 
a= aA Zo(f) (2.80) 


where tan 6 is the dielectric loss tangent and f is in GHz. 
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Table 16. MULTI-SECTION MICROSTRIP THEORETICAL RESULTS 


Ree Gem See Ne 
Static Effective Permittivity , 
Frequency-Dependent Eff. Permittivity ef, ( f) (2.6) 
Frequency-Dependent Impedance Z( f) 
Fide em) Semen @ 
Static Effective Permittivity ey; 
Frequency-Dependent Eff. Permittivity ¢),(f) 4.100 

Shunt Capacitance oo 
[es em seamen) 
Static Effective Permittivity ah 
F-equency-Dependent Eff. Permittivity a AAO) 

Shunt Capacitance oo 


Step Discontinuity Circuit Parameters 


Shunt Capacitance c?) Sa —— pF 
Shunt Capacitance c®) 0.213 pF 
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The total propagation loss (in dB) for each microstrip line is modeled in the circuit 


model by a resistive attenuator 7-networks. Figure 6.20 shows the equivalent circuit 


Table 17. MULTI-SECTION MICROSTRIP THEORETICAL LOSSES 


Narrow (0.8 em) Segment No.@) [| 


Dielectric Loss al) (2.10) 0.450 dB/A, 
ayy | coma 
Total Segment Loss ater 5 ee FT 0.487 dB/A, 
Wide (Sem) Segment No @) | | 
Dielectric Loss a (2.10) 0.484 dB/A, 
ay | on a, 
Total Segment Loss [ol Segment toss | «| ss aay 0.503 dB/A, 
[Wide (sem) Segment Re @ | 
am puna, 
Conductor Loss “ (2.8b) 0.021 dB/A, 


Total Segment Loss al?) Ed 0.498 dB/A, 


model for the microstrip filter of Figure 6.19. 















b. Simulation Results 
The microstrip filter was fabricated on G-10 epoxy dielectric material. 
Scattering parameter measurements were taken using the network analyzer. A 
piecewise-linear approximation of the measured excitation signal is used in the 
PSPICE circuit simulation. The transient analysis of the equivalent circuit model 
provides both reflection and transmission responses. Shown in Figure 6.21 is a com- 
parison between the measured and simulated reflection responses. The amplitude of 


the simulated response indicates that the theoretical losses are slightly greater than 
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than the measured losses. The magnitude spectrum of the reflection responses is 


shown in Figure 6.22. 





Figure 6.20 Multi-section equivalent circuit model 


The multi-section microstrip filter is also simulated by PUFF, a computer- 
aided design (CAD) program for microwave integrated circuits [Ref. 55]. The PUFF 
simulation solves the scattering parameters for the microstrip filter described in Fig- 
ure 6.19. The inverse FFT is performed on the scattering parameters to obtain the 
reflected and transmitted impulse responses. The PUFF program input excitation 
is a gaussian-shaped unit amplitude pulse. In order to accurately compare both the 
simulation results, a piecewise-linear approximation of the PUFF excitation was 
used in the PSPICE transient analysis. Figures 6.23 and 6.24 compare the reflec- 
tion and transmission responses, respectively. The PUF:* simulation model does 


not account for propagation losses. 














Figure 6.24 illustrates the effect of losses in the equivalent circuit model. The 
influe ..e of loss and dispersion on the multiple reflected signals can be observed in 


the Figure 6.23 by the apparent phase shift in the PSPICE simulation. 
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Figure 6.21 Comparison of PSPICE reflection response waveforms 
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Figure 6.22 PSPICE reflection response magnitude spectrum 
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Figure 6.23 Comparison of simulated reflection responses 
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2. Layer-Probing Simulation Results 
A lumped-distributed equivalent circuit is developed for a three-sectioned 


microstrip step discontinuity structure. A PSPICE transient analysis is performed 
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Figure 6.24 Simulated transmission responses 


using the equivalent circuit model. Reflection and transmission response data will be 
applied to the layer reflection coefficient estimation and layer-probing algorithms. 
The initial reflection impulse response data is applied to a Fortran program im- 
plementation of the layer reflection coefficient estimation algorithm. As shown in 
Figure 6.25, a sign change in the layer reflection coefficient sequence corresponds 


to an abrupt width change on the microstrip line. An estimated group velocity of 
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2.53 mm/sample is computed from a measured propagation time delay of 707.9 ps 


and a structure length of 112 mm. Figure 6.26 shows the impedance profile 


2 
a 9 2 92 908 G08 20-0 6-08 00-8 G0: 
é 


_ 
= 
Ld 
3 
ve 
Lo 
uJ 
oO 
O 
z 
° 
= 
O 
uJ 
J 
L 
LJ 
8 
od 
uJ 
5 


20 40 60 
AVERAGE LENGTH ( 2.533 mm/sample ) 





Figure 6.25 Multi-section layer reflection coefficient sequence 


computed from the estimated layer reflection coefficients. The accuracy of the 
impedance profile degrades rapidly beyond the first discontinuty because of loss and 
dispersion. Nevertheless, we are able to detect the relative changes in the impedance 
profile and then estimate the number of discontinuities within the structure. This 


result is more qualitative than quantitative. The detection resolution performance 
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of the algorithm is dependent upon the sampling interval used. Accurate layer 
reflection coefficients are obtained only when the sampling interval is sufficiently 
small. 
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Figure 6.26 Multi-section impedance profile 


The WLS algorithm was used to efficiently compute a 1st-order AR model 
from the reflected impulse response shown in Figure 6.27. The layer-probing algo- 
rithm next formulates an appropriate inverse filter using the estimated AR model 
parameter. A new excitation signal is simulated from the difference equation de- 
scription of the first inverse filter transfer function. We create a piecewise-linear 


approximation of the customized input and initiate a second transient analysis in 
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PSPICE to simulate the pulse propagating in the microstrip. The resulting second 


reflection response is windowed to eliminate the reflected energy caused by the first 
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Figure 6.27 Initial reflection response waveform 


discontinuity. The windowed reflection response is shown in Figure 6.28. The reflec- 
tion response is now deconvolved with the estimated transmission impulse response 
of the first discontinuity. Figure 6.29 shows the result of the deconvolution. The 
WLS algorithm is used to estimate another 1st-order AR model that will describe 
the second discontinuity. The previous AR model parameter estimates are used to 


design a 2nd-order inverse filter. Next, we excite the difference equation description 
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of the inverse filter using the original input pulse. A third excitation sequence is gen- 
erated to probe the microstrip. Once again a piecewise-linear approximation of the 


third customized input is entered in PSPICE. Figure 6.30 shows the transmission 
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Figure 6.28 2nd windowed reflection response 


response simulated by the multi-section equivalent circuit model. The simulated 
transmission response immediately reveals the transmitted impulse response of the 
third discontinuity. The transmitted impulse response is identified as the exponen- 


tially decaying transient shown in Figure 6.30. Finally, a third AR model is 





110 








obtained from the transmitted impulse response using the WLS algorithm. The 
purpose of the layer-probing algorithm was to resolve the reflected and transmitted 


impulse response for each cascaded step discontinuity. This was accomplished 
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Figure 6.29 2nd reflected impulse response after deconvolution 


by iterately “probing” the multi-section structure with customized signals synthet- 
ically produced by inverse filtering. Table 18 compares the estimated AR model 


parameters with the theoretical values used in the simulation. 
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The multi-section microstrip structure of Figure 6.19 was excited by 


a 18.5 GHz incident pulse. The performance of the layer reflection coefficient 


Table 18. MULTI-SECTION MODEL PARAMETER ESTIMATES 
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Figure 6.30 Third discontinuity transmitted impulse response 
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estimation algorithm is studied using the measured reflection response shown in 


Figure 6.31. The layer reflection coefficients are shown in Figure 6.32. 
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Figure 6.31 Measured reflection response due to the 18.5 GHz pulse 


However, the estimated reflection coefficient values rapidly decreased in amplitude 
after the 30th sample value. In Figure 6.31, we can observe a significant initial 
reflection which is followed by several lower amplitude transients. The primary 
reflection is due to the first discontinuity, while the subsequent transients are dis- 
torted by the multiple reflections. Next, the impedance profile is computed from 


the estimated reflection coefficients using an initial impedance value of 502. 
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Referring to Figure 6.33, a characteristic impedance of 12 is estimated for the first 


step discontinuity as compared with a theoretical value of 9.50 at 4 GHz. 
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Figure 6.32 Layer reflection coefficient estimate using measured data 


The impedance profile estimates a 50M microstrip line following the initial disconti- 
nuity. This narrow impedance estimate compares favorably to the theoretical value 
of 532. 

The simulation and measured results presented in this section have demon- 
strated the performance of the layer-probing algorithm. However, the algorithm 
assumes that the multiple discontinuities can be accurately modeled by cascading 


minimum-phase transfer functions. A practical realization of the algorithm could 
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be achieved by using optoelectronic picosecond sampling techniques and real-time 


digital signal processing. Finally, the layer impedance profile, using measured re- 
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Figure 6.33 Characteristic impedance profile for multiple discontinuities 


flection data, produced excellent impedance estimates for both the initial cascaded 


microstrip step discontinuity and the narrow 0.3 cm microstrip line. 
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Summarizing, we have presented both experimental and simulation re- 
sults for three microstrip test structures. First, pulse propagation simulations were 
performed using both the equivalent circuit and ARMA/AR parametric models of 
an IC microstrip interconnection. Time-domain measurements were compared with 
the simulation results. Second, a cascaded microstrip step discontinuity was inves- 
tigated. Using the equivalent circuit and parametric models derived in Chapter IV, 
reflection and transmission response simulations were compared with experimental 
results. Third, an equivalent circuit model was developed for a lossy multi-section 
microstrip structure consisting of several cascaded step discontinuities. The circuit 
model was validated by comparing PSPICE and PUFF simulations with microwave 
measurements. Finally, the parameter estimation algorithms, presented in Chapter 
V, were applied to simulated and measured response data in order evaluate their 


performance. 
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VII. SUMMARY AND CONCLUSION 


With the development of automated optoelectronic measurement systems suit- 
able for the manufacturing environment comes the problem of how to characterize 
microwave integrated circuit (MIC) devices and interconnections in MMIC’s. This 
thesis has presented new computer-aided design models and parameter estimation 
techniques to model these MIC interconnections using picosecond time-domain mea- 


surements. 


A. SUMMARY OF SIGNIFICANT RESULTS 

Equivalent lumped-distributed circuit models were derived for lossy dispersive 
microstrip transmission line structures. The circuit models consisted of resistive and 
reactive lumped-elements, and distributed lossless transmission-line sections. The 
models were compatible with general computer-aided design techniques and simu- 
lation programs, such as PSPICE. In addition, the circuit models for the cascaded 
microstrip step discontinuity accurately characterized the reflected and transmitted 
wave propagation in the structure. These models mark a significant improvement 
over empirical modeling for which accurate generalized closed-form solutions are 
not available for CAD programs. 

The picosecond sampling performance of current optoelectronic measurement 
systems motivates a digital signal processing framework for the modeling and anal- 
ysis of high-speed transient signals. Parametric models of lossy dispersive microstrip 
transmission-lines provide an innovative modeling approach not previously pre- 
sented by microwave engineers. The parametric model approach is emphasized 
for several reasons. First, the fitting of an autoregressive moJel to an observable 
time series is a linear process which can be handled using well tried and highly 
efficient computing algorithms. Secondly, an autoregressive moving-average model 


provides the smallest number of parameters for an optimal and most parsimonious 











representation of the propagation process. Finally, the identification problem in 
parametric modeling, in other words choosing an AR or ARMA representation for 
a microstrip line, was shown to be intimately related to the critical physical prop- 
erties. We have demonstrated how the effective microstrip permittivity of a MIC 
interconnection can be evaluated from ARMA or AR model parameters. Further- 
more, parametric model simulations accurately characterized the picosecond pulse 
propagation on lossy dispersive microstrip structures. A first-order AR model of 
the cascaded microstrip step discontinuity formed the basis for the development of 
the multi-section microstrip structure model. 

Several parameter estimation techniques were presented. Both deterministic 
and stochastic algorithms were developed to estimate the ARMA and AR model 
parameters. The weighted least squares (WLS) algorithm estimated the ARMA dig- 
ital filter coefficients or model parameters from a prescribed unit impulse response. 
The algorithm rapidly converged when applied to several finite-length deterministic 
examples. 

A discrete transformation of the network function approximation method, in- 
troduced by Elmore [Ref. 46], was presented. This method approximated either a 
first-order “all pole” or a second-order “pole-zero” network function from measured 
transient pulse delay and rise times. A first-order AR model and a second-order 
ARMA model were derived using the impulse invariant transformation of the re- 
spective network functions. 

A layer reflection coefficient estimation algorithm was derived. This algorithm 
estimated the reflection coefficients from the sampled reflection response of a lossy 
multi-section microstrip structure. The algorithm assumed a constant attenuation 
factor and group velocity throughout the structure. The algorithm detected major 
discontinuities within the structure. 

Finally, the layer-probing algorithm was developed to characterize multiple dis- 


continuities within a multi-section microstrip transmission-line, consisting of several 
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step discontinuities. A minimum-phase parametric model described each discontinu- 
ity. Inverse filtering and deconvolution techniques were employed by the algorithm 


to iteratively probe successive sections of the structure. 


B. CONCLUSIONS 

The research ccnducted for this thesis has shown that the parametric modeling 
approach provides an alternative framework to characterize and analyze dispersive 
pulse propagation on microwave integrated circuit interconnections. The model- 
ing of dispersive propagation on silicon and gallium arsenide microstrip structures 
is important because it points to a limitation in IC speeds even when transistor 
switching speed continue to improve. The measurements conducted and reported 
in this work have verified that the proposed equivalent circuit models are compati- 
ble with available computer-aided design (CAD) programs. The circuit models were 
developed using a combination of frequency-dependent expressions available in the 
literature and original analysis by the author. 

The development of optoelectronic measurement systems will provide the ca- 
pability to directly sample high-speed transients on the semiconductor substrates 
without introducing significant parasitics. With the availability of time-domain 
measurements, we developed new signal processing algorithms to estimate both the 
ARMA and AR model parameters. 

The layer-probing algorithm presented a technique to iteratively identify indi- 
vidual step discontinuities in a multi-section transmission-line. Simulation results 
showed how destructive interference problems caused by dispersion and multiple 
reflections degraded the algorithm for some structures investigated. This problem 


was overcome by measuring and modeling the transmitted impulse response. 
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C. FUTURE DIRECTIONS FOR RESEARCH 

The objectives of the research were successfully accomplished. Several sugges- 
tions for future research are now presented. 

Researchers are currently studying optical techniques to measure nonlinear and 
multiport characteristics of MMIC components and circuits. Based on these experi- 
ments, nonlinear parametric models should also be investigated. The high frequency 
sampling capability of photoconductor devices can be integrated with pulsed laser 
technology to develop a picosecond time-domain measurement system using a white 
noise excitation. Here, several stochastic-based parameter estimation algorithms 
can be applied to the modeling problem. 

Equivalent circuit and parametric models need to be developed for coupled 
microstrip, coplanar waveguide and other microstrip discontinuities, such as bends, 
crossings and the asymmetrical step discontinuities. 

Improvements to the layer-probing algorithm should include: 1) the develop- 
ment of a method to estimate the attenuation factor of each lossy section. 2) the 
implementation of a minimum-phase ARMA parametric model for each section, and 
3) the integration of cross-correlation techniques. 

Another important aspect of this research is its extensions to other engineering 
problems. Future efforts envision the application of the proposed system identifica- 
tion algorithms to ultra-wideband (UWB) radar target identification. Emphasis is 
required on time-domain electromagnetics in measuring and computing the impulse 
response of targets, radar-absorbing materials, transmitting and receiving antennas. 


and various other media. New signal processing algorithms must be developed for 


impulsive radar. 





Appendix A_ Derivation of the Schur Recursion Algorithm 


A recursive in order solution for the AR lattice filter reflection coefficients is 


presented. This derivation of the so-called Schur algorithm will show the computa- 


tional efficiency of this algorithm over the Levinson-Durbin algorithm [Ref. 12: pp. 


6-22]. The algorithm derivation will assume that the signal data is stationary. 


1. List of Symbols and Definitions. 


The following symbols are used in the derivation: 


y(n) 


The autoregressive (AR) signal sequence. 

The prediction error filter coefficient vector. 

The forward prediction error. 

The backward prediction error. 

Expectation operator. 

The autocorrelation matrix of the signal y(7). 

The z-transform of the forward model parameter. 
The z-transform of the backward model parameter. 
The (p+ 1)-th order model reflection coefficient, or lattice filter 
reflection coefficient. 

The z-transform of the forward prediction error. 
The z-transform of the backward prediction error. 


The (p)-th order normalized forward prediction error. 


(p)- 
The (p)-th order normalized backward prediction error. 
he (p)-th order Schur forward recursion parameter. 


T 
The (p)-th order Schur backward recursion parameter. 








2. Algorithm Derivation. 
The present output data y(n) can be estimated from the weighted summation 
of the past values of the output for the (p)-th order autoregressive (AR) model. This 
estimate is expressed as [Ref. 46: p. 150] 


Pp 
y(n) = S> by y(n —2) (A.la) 
i=} 
or in matrix form 
y(n) = b' y(n) (4.16) 


where 


y? (n) = [y(n — 1) y(n — 2) ... y(n — p)] 


bT = [by by... byl 


The prediction error is defined as the difference in the signal and its estimate or 


= y(n) — 4" y(n) (4.2) 


The mean squared prediction error, expressed as a function of the prediction error 


filter weights 6, is formulated as (Ref. 56: p. 20] 


E{e(n)} = 6 [Ryy] b-— 207 Ry, + Ryy(0) (A.3) 
where 
Ryy(0) Ryy(-1) eee eee Ryy(—p + 1) 
Ryy() Ryy(0) «..-. Ryy(p) 
[Ryy] = . iz . 
Ryy(p = 1) Ryy(1) Ryy(0) 
and 


R= yl) Ry) ce Rel) 


The solution of (A.3) for the optimum filter weights is given by 


Bop = [Ryy) Ryy (As) 


which is commonly call-? the normal equations [Ref. 56: p. 22]. The normal 


equations can be written in matrix form as 


Ryy(0) Ryy(—1) ee ecee Ryy(—p + 1) by Ryy(1) 
Ryy() Ryy(0) «ee. Ryy(p) bo Ryy(2) : 
: : om : : = : (4.5) 
Ryy(P ==1) oes Ryy(1) Ryy(0) bp Ryy(P) 


The autocorrelation matrix [R,,} is Toeplitz if the signals are stationary which is 
assumed in this derivation. The Levinson algorithm provides the solution of the 
(p+1)-th order model parameters by using the previously determined p-th order 
model parameters. [Ref. 46: pp. 81-85 ]. The (p+ 1)-th order model parameters 


are obtained recursively by 


BlP+})(z) = B'?)(z) - z—1 pl? +1) BY) (2) (A.6a) 
BOHN) 2) = 2-1 BW) (z) — k(P+1) BlP) (2) (A.6b) 


where B'?)(z) is the p-th order forward model parameter transfer function, B‘)(z) 
is the p-th order backward model parameter transfer function, and k{?+!) is called 
the reflection coefficient. The z-transform of the forward and backward prediction 


error is obtained by multiplying (A.6) by Y(z). 


EPA (2) = Bl) 2) — 271 RP th) FP) 2) (A.7a) 
EPM (2) = 21 B)(z) — Ret) BlP)( 2) (A.7b) 


In the time domain the prediction error equations become 


e(P+1)(n) .f e'P)(n) -m klP+ 1 ElP) (n —1) (A.8a) 
BPH) (n) = @)(n — 1) — k(P+1) 6?) (2) (A.8b) 


To find e'?t)(n) and é?t+))(n) from e'?)(n) and é'?)(n), respectively, it is only 
necessary to determine k'?+!), Since the zero order prediction of a signal is y() = 0, 
the zero order error is 


e°(n) = @(n) = y(n) (A.9) 











and equations (A.8a) and (A.8b) can be realized as a lattice structure as shown in 


Figure A.1. 
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Figure A.1 Second Order AR Lattice Filter Model 


The reflection coefficient k{?+!) can be determined by minimizing the square of the 
expected value of equation (A.8a) and (A.8b) with respect to the forward (ky) and 
backward (ky) reflection coefficients, respectively. The forward reflection coefficient 
is 


EfelP)(n)e)(n a 1) 


(PN Ne (A.10) 
E {a(n — 1)8)(n — 1)} 
and the backward reflection coefficient is 
Ed el) (ne) (n — 1) 
kyPt) = E{enjen - 1} } (A.11) 


E{ a°)(n)e()(n)} 


The forward and backward reflection coefficients can be related to the geometric 
mean form as (Ref. 57: pp. 434-437] 


E{e'?)(n)é(P)(n — 1)} 


VEX [e(n)]? JEL [En —1)]?} 


k(P+}) — (A.12) 
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By substituting (A.12) into (A.10), the forward reflection coefficient becomes 


E{ [el)(n)]” JE{ [Pn — 1)]7} 


fees ph) Ne ee A.13a 
, E{ [e\(n -1)}"} aoe 
eyPH)) = girs) __ lend (A.13b) 
I| [2(n — 1)" |] 
where 
lle™(n) |] = Ef [e(n)]” } 
and 


Zn — 1) = VEL 2(n - 1)]} 


Similarly, the backward reflection coefficient (A.11) yields 


kt) = gira VEEP eC len — WT} (A.14a) 
E{ [elP)(n)]"} 
kyPt)) = porn lem —DIl (A.148) 
[en] I | 


Substituting (A.13b) and (A.14b) into (A.9a) and (A.9b), respectively yields 


(p) 
elPt+1)(n) = e!?)(n) = K(P+D el?) (n es je ete (A.15a) 
Il [2m —1)]" | 
(P+) (K) = Bn — 1) — PHP (y oe (4.158) 
I [e(n)} "I 
The normalized (p + 1)-th order forward and backward — errors are defined 
as (p41) 
PHN (ny) = A) 
é (n) = [e(n)i (A.16a) 
é(P+1)(m — 1) 
pe ap ak A.16b 
r (n jz7(n— DI ( ) 
Factoring ||e?)(n)|| from the right side of (A.15a) and normalizing yields 
(Pp) 
=(p+1) ae lle‘?(n)| =(p) (pt )a(P)(, 7 
kee era Es " i} aie) 
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Factoring |é?)(n — 1)|| from the right side of (A.15b) and normalizing yields 


AP) (n — 
ee ae (cies ass 81 eT 7 
el are cer {i (n — 1) — Rete (n)} (A.17b) 
We will next show that 
lle+1)(n)]] 2 
ih [A(e+4)] 
= \2(n — 
eP)(n — 1 1 
ial aa (A.18b) 
(p+i)(, — 
eee ON ha. [korn]? 
Proof: Equation (A.18) is defined as 
1 
2 vi] 
(p) 
Un) le { {een 
eae Ski een aries: Lee (A.19a) 


feP*D(n)IP 


lef [eP*D(n — yy)? \)' 


Squaring (A.9a) and taking the expectation yields 
2 2 
e{ [e?*(n)] \ = e{ [e'")( (n)| Ve akPHEf ein yal”(n = iy} 


+ [RetDy €{ [ein - n) \ (A.196) 
Substituting (A.10) for ky+!), and simplifing 


e{ [e*)(n)] \ = e{ [en] } = bp PML elP)(nyeP)(n 2 1)} (A.19¢) 


Returning to equation (A.19a), we substitute (A.19c) to yield 


len) 1 
[ern] T 
e{ [et(n)}” \ - kp PAE l)(n)e)(n = 1)} 2 
= (A.19d) 


i 
kg PDEL eP)(n)e)(n = 1)} 2 


e{ [e((n)]? \ 


= 














Finally, we substitute (A. 13a) for kp) into (A.19d), and multiply the denomina- 
tor term e{ le (P) (n)}- \ by the ratio 


e{ [e(n)]” [a(n = 1)| : \ 
e{ [en] [a(n = ») \ 


simplifying yields 


Me IL 1 | 
ePFN(n) I~ ge Aa 
. 3 2 
PEL elP)(nye(n = )} 
e{ [e)(n)]? [een —1)]? } 
el ¢ (n 
Hen)” aie 


le(P+)) (nil 
FemMEH AY 


Following the same approach using equations (A.9b) and (A.11), it can be shown 
that 


alP)(m — J 
eer = TT a — = (A.20) 
i [aren 
The normalized lattice recursion equations are from (A.17a) and (A.17b) 
aad = 1 1 — k(p+1) é(P)(n) 
ae _ i [k(p+1)}? — (p+) 1 F(P)(n - 1) (A.21a) 
or aoe 
e+) (n) | _ prp(ptiy) | @?)(n) 
ea = [Kir rim — 1) (A.216) 
where 
1 (p+1) 
.(p+1)]) — 1 —k 
sg afi [k(p+1)}2 | _gben l | 








In order to obtain the Schur recursion equations we will introduce two new param- 


eters, a?)(n) and Bn —1) as 





a)(n) = E{2(k) y(n)} (A.22a) 
e(?)(n) y(n) 
=e{-_ TOIT — (A.226) 
Ryy(n) 
a®)(n) = Buy 226 
oe) = Ty) ore 


where R,,(n) is a (p x 1) correlation vector, and the norm |\y(n)|| = \/Ryy(0). 


Similarly, 
B(k)®(n — 1) = E{F)(n ~1) y(n )} (A.23a) 
eens (A208 

B)(n —1) = aT (A.23c) 


Using the normalized lattice recursion equation (A.21) and the previously defined 


Schur recursion parameters (a'?)(n)) and (Bp? Yin — 1)), the Schur recursion equation 


| Ny) > are | al?)(n) | (A.24) 


1S 
BPH (n B")(n — 1) 


We will next show the ((p+1))-th order reflection coefficient in terms of the 


ratio of Schur recursion parameters. From equation (A.22a) and (A.23a) we form 


an) E{e™(n) y(n) 
B(n-1) Ef (e(n —1) y(n)} 
__ (10 
~ ef rtel(n —1)} 


(A.25a) 


(A.25b) 


128 








Ese) (n) SNe (n —1 
_ Efe*M(n) hen — 1) nae 


E{z0P)(n — 1) ben) 
a a een en-f (A.25d) 
BP'(n -1) €{ [el(n)]? by /e{ [en -1)]”} 


Using stationarity and é{ [e”(n - 1») } = e{ (en) \, we have [Ref. 46: p. 
168] 


a(n) _ Efe(n) e(n —1)} 


TT A.25 
B®)(n —1) e{ [e(n - 1)]"} ( f) 
= E(P+)). (A.25g) 


3. Normalized Lattice Filter Realization. 
This section describes the normalized lattice filter as a direct method of 
solving k(t!) given the autocorrelation sequence (R,y,). Assume the zero order 


Schur parameters are 


a1) = R(1) 
B)(0) = R(0) 


Inserting these initial values in (A.25) we have 


,) = a)(1) 
B()(0) 

_ Ril) 

= Fi) 





The first order Schur parameters for (n = 1,2) are solved using the recursion equa- 
tion (A.24) as 


aaa [a!(1) = EO ((0)) 


a1) = 














and 


a()(2) = vende (0) (2) — (1) 8)(1)] 
Oy 


where a!)(2) = R(2) and 6((1) = R(1). The B“)(n) parameters are 


Wits : Kay? 


1 


Vi - [ROY 


where a!°)(2) = R(2) and 8(0)(1) = R(1). Substituting the previous results into 
(A.25g) yields 


6% (0) = = [A°(0) - — kM g(0)/ (1)] 
and 


B(1) = (ac) 7 KMa((2)| 


42) — al)(2) 
B)(0) 

_ R(2)—kO RQ) 

~ RO) —kO ROA) 


The Schur recursion equation can be realized as a p-th order normalized lattice 


structure as shown in Figure A.2. 
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Figure A.2(a) Zero Order (n = 1) Normalized Lattice Filter 
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Figure A.2(b) First Order (n = 2) Normalized Lattice Filter 
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APPENDIX B_ PSPICE EQUIVALENT CIRCUIT LISTINGS 


IC INTERCONNECTION (PI-NETWORK AND MATCHED IMPED TERMINATION) 


* MICROSTRIP TRANSMISSION LINE SIMULATION (LCL) Pl - MODEL 
* TEST STRUCTURE ( LENGTH « 500E-06 m) 

* SILICON SUBSTRATE ( ER = 11.7) 

.OPT ACCT LIST NODE OPTS NOPAGE RELTOL«.001 
.WIDTH IN#2133 OUT=80 

OPTIONS LIMPTS =» 20000 

-OPTIONS ITLS = 50000 

OPTIONS TRTOL:5 

.OPTIONS ITL4=25° 

V 10 PWL(0 0 3.75p .11 7.5p .25 11.25p 1.45 15p 1.9 18.75p 2 22.5p 1.9 26.25p 1.85 37.5p 1.6 52.5p 1.2 
75p .72 78.75p 0) 

RS 12 74.23 

T1203 0 2074.23 TD=1.25ps 

R12 0 1.065K 

R223 10.4 

R33 0 1.065K 

C130 0.03pF 

L1340.332nH 

C2 40 0.03pF 

RL 40 74.23 

.TRAN 1.25pS 80ps 0 1.25pS 

.PROBE 

.END 


IC INTERCONNECTION (T-NETWORK AND MATCHED IMPED TERMINATION) 


.OPT ACCT LIST NODE OPTS NOPAGE RELTOL+x.001 
-WIOTH IN#133 OUT=80 

OPTIONS LIMPTS # 20000 

OPTIONS [TLS = 20000 

OPTIONS TRTOLe5S 

OPTIONS [11405 

* PIECEWISE-LINEAR APPROXIMATION OF NORMALIZED PCE SAMPLER 1 
V 10 PWL(0 0 3.75p .11 7.5p .25 11.25p 1.45 15p 1.9 18.75p 2 22.5p 1.9 26.25p 1.85 37.5p 1.6 §2.5p 1.2 
75p .72 78.75p 0) 

RS 1 2 74.23 

R120 1.065K 

R223 10.4 

R330 1.065K 

T1203 0 Z0#74.23 TDe2.23ps 

C130 0.06pF . 

T230 40 20=874.23 TD=2.23ps 

RL 40 74.23 

.TRAN 1.25pS 100pS 0 0.05ps 

PROBE 

.PRINT TRAN V(4) 

.END 
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CASCADED-STEP DISCONTINUITY PSPICE PI-MODEL SIMULATION 


* DIELECTRIC: (G-10 ER = 4.2) 


* INPUT SIGNAL: (PIECEWISE-LINEAR APPROXIMATION OF THE MEASURED 
* GAUSSIAN PULSE OF THE HP8510 NETWORK ANALYZER) 


OPT ACCT LIST NODE OPTS NOPAGE RELTOL=.001 


-WIDTH IN=133 OUT=80 
-OPTIONS LIMPTS = 20000 
-OPTIONS ITLS5 = 50000 
-OPTIONS TRTOL=5 
OPTIONS [T1425 


V 10 PWL(0 0 100P .06 300P .17 400P .37 800P 1.88 900P 2.02 1NS 1.93 1.4NS .41 1.5NS .17 1.6NS 


044 1.7NS 0) 

RS 12.41.24 

T1 2030 2004124 TD=750PS 
R1303.118K 

R23 41.091 

R3 40 3.118K 

T2 405.0 20041.24 TD=214PS 
C1507.13PF 

L4 5 6 0.362NH 

C26 07.13PF 

T3 6 0 70 20041.24 TD=214PS 
R4703.118K 

R578 1.091 

R66 03.118K 

RL 8 041.24 

TRAN 100PS SNS 0 30PS 
PROBE 

.PRINT TRAN V(2) 

.PRINT TRAN V(8) 

.END 
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CASCADED-STEP DISCONTINUITY PSPICE SHUNT CAPACITOR MODEL SIMULATION 


* DIELECTRIC: ( EPOXY G-10 ER = 4.2) 
* INPUT SIGNAL: (PIECEWISE-LINEAR APPROXIMATION OF THE MEASURED 
* GAUSSIAN HP8S10 NETWORK ANALYZER SOURCE) 


OPT ACCT LIST NODE OPTS NOPAGE RELTOL«.001 
WIDTH IN=133 OUT=80 
OPTIONS LIMPTS = 20000 
OPTIONS ITLS5 = 50000 
OPTIONS TRTOL=S 

OPTIONS [T1425 

V 1 0 PWL(0 0 100P .06 300P .17 400P .37 800P 1.88 900P 2.02 1NS 1.93 1.4NS .41 1.5NS .17 1.6NS 
044 1.7NS 0) 

RS 12 41.24 

T1203 0 20=41.24 TDx750PS 
R13 02316K 

R23 41.409 

R3 40 2.316K 

T2 405 0 20041.24 TD#214PS 
C15 0 10.91PF 

T45 0 6 0 20m7.424 TD=24.38PS 
TS 6070 20=7.424 TD=24.38PS 
C270 10.91PF 

T3 7080 20041.24 TD=214PS 
R48 02.316K 

R58 9 1.409 

R6 9 0 2.316K 

RL 90 41.24 

-TRAN 100PS SNS 0 30PS 
-PROBE 

-PRINT TRAN V(2) 

PRINT TRAN V(9) 

-ENOD 
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MULTI-SECTION MICROSTRIP FILTER LAYER-PROBING PSPICE SIMULATION 


* DIELECTRIC:EPOXY (ER = 4.3), H=1.59mm, DESIGN FREQ = 8.79GHz 
* CHARACTERISTIC IMPEDANCE = 59.2 Ohms (THEORETICAL) 

OPT ACCT LIST NODE OPTS NOPAGE RELTOL=.001 

-WIOTH IN=133 OUT=80 

-OPTIONS LIMPTS = 100000 

OPTIONS ITLS = 100000 

OPTIONS TRTOL«5 

OPTIONS [T1405 





* THE FOLLOWING PIECEWISE-LINEAR APPROXIMATION DESCRIBES AN INITIAL 100PS 
* GAUSSIAN PULSE HAVING A 3DB BANDWIDTH OF 8.79 GHZ 
“V3 10 PWL(0 0 10P .03 20P .33 30P .968 40P 1.68 SOP 2 60P 1.68 70P .968 80P .33 90P .03 100P 0) 


* THE FOLLOWING PIECEWISE-LINEAR APPROXIMATION DESCRIBES THE SIMULATED 
* INVERSE FILTER OUTPUT USING THE 1ST DISCONTINUITY AR COEFFICIENT. 

°V3 10 PWL(0 0 10P .03 20P .302 30P .668 40P .8 50P .47 60P -.1426 70P -.562 BOP -.552 SOP -.274 
100P 0) 

* THE FOLLOWING PIECEWISE-LINEAR APPROXIMATION DESCRIBES THE SIMULATED 
* 2ND-ORDER INVERSE FILTER OUTPUT USING 1&2 DISCONTINUITY AR COEFFICIENTS. 
“V3 10 PWL(O 0 10P .03 20P .276 30P .398 40P .206 50P -.238 60P -.558 70P -.436 80P -.053 90P .22 
100P .212 110P .024 120P 0) 

RS 1259.22 

T0 2030 Z0e59.2 TD=1P 

R130 831.7 

R23.48.48 

R3 40 831.7 

T2 405 0 Z0=59.22 TD=292.6PS 

C15 03.57PF 

T4 5 0 6 0 Z0859.22 TDa77.7PS 

R460 1.55K 

RS 674.52 

R6 7 0 1.55K 

TS 70 8 0 Z0=59.22 TDa77.7PS 

C28 0 1.5PF 

T6 80 9 0 Z0=59.22 TD=77.7PS 

R790 1.55K 

R8 9 10 4.24 

R9 10 0 1.55K 

T7 10 0 11 0 Z0#59.22 TD=77.7PS 

C3 1103.57PF 

T8 11 0 12 0 Z0u59.22 TD=104,4PS 

R13 12 0 2.323K 

Ri4 12 13 3.02 

R15 13 0 2.323K 

RL 13 0 59.22 

TRAN 5PS 2.5NS 0 5PS 

PROBE 

PRINT TRAN V(2) 

.PRINT TRAN V(13) 

.END 
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